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ABSTRACT  ' 

Whitham’8  results  for  slowly  varying  wavatrains  ara 
ganaralizad  to  include  (i)  explicitly  the  effect  of  modal 
dependence t (ii)  the  effect  of  low-order  linear  or  honlinear 
nonconservative  terms  and  their  role  in  modifying  the  basic 
. competition  between  frequency  and  amplitude  dispersion, 

(iii)  the  effects  of  moving  media,  rotational  or  irrotational, 
and  space-time  inhomogeneities,  and  (iv)  the  effect  of 
high-order  dispersive  or  diffusive  modifications  to  the  low-order 
amplitude  and  phase  equations.  A "modal  law"  is  also  derived 
which,  when  integrated  over  the  cross-space,  leads  to  a 
generalized  Whitham  action  law  modified  by  a source  term  that 
is  introduced  by  nonconservative  effects i the  use  of  complex 


frequencies  here  is  avoided,  thus  enabling  linear  and  nonlinear 

problems  to  be  treated  on  an  equal  basis.  A self-consistent 

\ 

perturbation  method  is  presented  that  leads  to  generalizations 
and  extensions  of  those  ideas  due  to  Whitham;  Hayes,  Davey, 
Stewartson,  Bretherton  and  Garrett,  Taylor,  Landahl,  and 
others,  and  shows  how  each  of  these  generalized  subsets  appears 
as  a specific  limit  of  a broader  theory.  A number  of  examples 


Illustrating  the  basic  ideas  are  presented  and  which  lead  to 


simple  formulas  that  can  be  used  in  many  direct  applications 
of  the  general  theory.  Some  qualitative  ideas  on  the  basic 
effect  of  the  derived  high-order  terms  are  given,  especially 
with  regard  to  the  question  of  kinematic  shocks,  and  their 
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basic  rslationshlps  with  other  hlgh-order  theories,  for  example, 
Stewartson  and  Stuart’s  parabolic  diffusion  equation  for  wave 
amplitude,  and  a particular  nonlinear  Schroedinger. equation 
that  arises  In  studies  of  hlgh-order  wave  dispersion,  are  also 
discussed. 

Based  on  the  above  results,  some  new  Ideas  on  wave 
stability  are  presented.  One  particularly  Important  example 
discusses  Landahl's  focusing,  trapping  and  amplification 
mechanism  for  waves  riding  on  Inhomogeneous  flows  In  light  of 
an  analogy  with  the  Instability  of  transonic  nozzle  flows,  and 
some  speculative  Ideas  on  kinematic  shock  formation  are  given. 

The  discussion  essentially  parallels  the  simple  arguments 
originally  put  forth  by  Kantrowitz  in  gas-dynamics  but  in  the 
wave-mechanical  context;  however,  detailed  mathematical  analysis 

of  the  entire  Instability  mechanism  Is.  not  pursued  In  the 

.1 

present  paper.  Pursued  In  some  detail  here,  however.  Is  the 
study  of  wave-mechsmical  discontinuities  (which  are  presumed 
to  arise  out  of  some  wave-llke  Instability),  both  with  and 
without  the  effect  of  mean  flow.  Some,  simple  calculations  are 
given  for  gravity  auid  capillary  waves  that  Illustrate  the 
general  approach  undertaken,  and  a general  theory  Is  also 
presented  that  deals  with  wave -mechanical  shock  stability  as 
would  be  Influenced  by  low-order  nonllpearltles  and  Inhomogeneltles 
in  the  mbdlum;  the  theory  essentially  sets  the  foundation  for 
a general  wave-mechanical  approach  that  closely  parallels 
analogies  developed  from  nonlinear  acoustics.  Other  examples 
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treated  include  simple  generalizations  of  instability  theories 
due  to  Whitham,  Taylor » and  others,  and  a newly  discovered  (and 
somewhat  speculative)  anti-cascading  phenomenon  arising  from 
the  idea  of  high-order  "negative  wave  diffusion". 


I.  Introduction. 


Perhaps  the  most  exciting  area  in  applied  mathematics 
research  today  is  in  the  field  of  nonlinear  wave  mechanics  and 
its  application  to  the  wave  stability  of  continuum  systems. 

The  greatest  recent  advances  can  be  found  in  the  mid«1960*s, 
as  we  shall  see,  and  many  of  the  initial  advances  and  applications 
were  spearheaded  by  V^Thitham  and  his  coworkers.  The  subject,  as 
it  stands  today,  is  founded  on  the  contributions  of  many 
investigators,  among  them  Hayes,  Stewart son,  Stuart,  Davey, 
Bretherton  and  Garrett,  Landahl,  Lighthill  and  others t the 
cumulative  body  of  knowledge  developed  by  these  authors 
is  capable  of  handling  low-order  frequency  and  amplitude 
dispersion,  modal  effects,  wave  damping,  high-order  wave  diffusion 
and  dispersion,  and  other  simultaneous  effects.  However,  it  was 
often  the  case  that  various  of  the  particular  theories  were 
devised  specifically  to  deal  with  particular  problems  and  their 
extension^ to  include  higher-order  effects  and  additional 
physical  features  could  not  be  made.  The  broad  generality  of 
an  idea  could  sometimes  not  be  readily  perceived,  this  situation 
leading  to  many  rather  specialized  and  sometimes  conflicting 
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theories.  The  basic  purpose  of  the  present  paper  is  to 
reconsider  anew  the  general  problem  of  describing  a slowly 
varying  wavetrain  by  constructing  a broader  unifying  theory 
that  bears  as  few  limiting  assumptions  as  the  physics  permits 
and  which  reduces,  in  specific  limits,  to  those  theories  referred 
to  above. 

Numerous  applications  of  more  or  less  "standard”  nonlinear 
techniques  (e.g, , multiple  scaling  methods,  averaging  methods, 
etc. ) have  been  made  in  recent  years  to  problems  involving  wave 
motions  in  continuous  media  and  an  extensive  survey  of  these 
studies  appears  in  Chin  (1976).  These  applications  are  mainly 
based  on  several  fundamental  papers  and  we  will  here  summarize 
and  discuss  the  basic  ideas  and  theories.  Probably  the  first 
application  of  WKB  averaging  in  fluid  mechanics  appears  in 
Whitham  (1962),  in  which  (for  finite  depth  gravity  waves)  the 
flu id -dynamical  conservation  laws  were  phase-averaged  so  as  to 
produce  the  modulation  equations  governing  coupled  large-scale 
changes  in  the  wave  and  mean  flow  variables.  In  this  averaging 
the  Ansatz  consisted  of  the  related  uniform  wave  solution,  and 
one  additional  postulate,  that  of  wave  conservation,  was 
introduced.  This  averaging  method  was  later  used  in 
Whitham  (1965a)  to  find  slowly  varying  wave  solutions  of  nonlinear 
equations  in  general  and  there  the  analogous  role  of  the  basic 
uniform  solution  as  compared  to  that  of  the  sinusoidal  wave  in 
linear  stationary  phase  theory  was  advanced.  Whitham  (1965b) 
in  a further  development  showed  how  the  same  modulation 
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equations  for  a system  of  given  equations  could  be  derived  by 
evaluating  the  related  Lagrangian  density  with  the  uniform 
wave  solution  (the  slowly  varying  solution  is  assumed  to  behave 
locally  like  the  planar  wave),  phase-averaging,  and  then  taking 
appropriate  variations  with  respect  to  the  wave  and  mean  flow 
parameters  (which,  during  the  averaging,  are  held  constant). 

The  power  of  the  method  is  seen  in  its  simplicity.  For  example. 
In  its  application  to  finite-depth  gravity  waves  (Whitham,  196?) 
the  averaging  reproduced  in  a simple  way  the  concepts  of 
•radiation  stress"  and  the  results  of  Longuet -Higgins  and 
Stewart  (1961)  obtained  on  the  basis  of  detailed  asymptotic 
analysis.  Whitham  also  showed  how,  for  sufficiently  deep 
water,  the  competition  between  amplitude  and  frequency  dispersion 

led  to  an  ill-posed  initial  value  problem  which  in  a sense 

1 

implied  instability!  the  modulation  equations  were  elliptic. 

For  shallow  water  the  governing  equations  were  hyperbolic,  and 
implied  a "splitting"  of  the  basic  wave  envelope.  These  ideas 
are  quite  general  and  apply  as  well  to  other  applications  of 
the  Whitham  theory.  We  add,  in  summarizing  this  paragraph, 
that ^Whitham* s results  and  his  "average  Lagrangian  method"  can 
be  derived  and  Justified  by  multiple-scaling  techniques,  for 
example,  as  discussed  in  Luke  ( !%(>  ) and  in  Whitham  (1970). 

..  While  Whitham* s variational  theory  was  successful  in 
describing  both  frequency  and  amplitude  dispersion,  it  was  not 
able  to  incorporate  the  effects  of  weakly  nonconservativs 
terms,  such  as  wave  damping.  This  failure  stemmed  from  the  lack 
of  a suitable  variational  formulation  and  thus  precluded  the 
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application  of  Whitham'a  ideas  to  many  practical  problems  (some 
recent  modifications  of  the  theory  (Jimenez  and  Whithara,  1976) 
now  permit  the  treatment  of  such  effects,  but  the  description 
is  confined  to  low-order  changes  and  does  not  appear  to  be  as 
general  as  the  theory  to  be  presented).  Moreover,  Whitham'a 
method  could  not  be  generalized  so  as  to  produce  the  high-order 
dispersive  terms  obtained  by  other  investigators  using 
direct  multiple-scaling  methods.  For  example,  a WKB  technique 
applied  to  the  study  of  slow  modulations  of  a Stokes  wavetrain 
showed  that  higher  order  terms  representing  modulation  rates 
and  not  nonlinearity  must  be  added  to  extend  the  validity  of 
Whitham's  equations  (Chu  and  Mei,  1970).  The  same  ideas  were 
reemphasized  by  Davey  (1972),  who  suggested  some  heuristic 
modifications  to  Whitham's  theory,  and  by  Davey  and 
Stewartson  (197^)*  who  considered  three-dimensional  wave 
packets  in  water  wave  applications.  The  apparent  discrepancies 
in  the  two  approaches  could  not  be  found,  as  noted  in  Hayes  (1973). 

For  example,  direct  multiple-scaling  methods  led  to  a nonlinear 
Schroedinger  equation  for  a certain  complex  wave  amplitude,  and 
this  high-order  equation  embodies  Whitham's  low-order  results 
(centered  wavenumber  expansions,  however,  are  required  in  its 
usual  derivation,  so  that  the  competing  effects  between 
amplitude  amd  frequency  dispersion  are  not  accounted  for). 

Furthermore  we  note  that  Whitham's  method,  as  it  presently 

) ^ I 

^ stands,  does  not  describe  high-order  diffusive  effects  of  the 

kind  described  in  Stewartson  and  Stuart  (1971).  which  are  j 

certainly  important  in  any  large-time  description.  Thus  one  j k 

has  in  addition  to  low-order  frequency  and  amplitude  dispersive 


terms*  those  terms  responsible  for  the  high-order  competition 
between  wave  dispersion  and  wave  diffusion,  (Recently 
Yuen  and  Lake  (1975)  showed  how  the  average  Lagrangian  method 
could  be  modified  so  as  to  produce  the  required  high-order 
dispersive  corrections,  however,  but  certain  implied  assumptions 
are  made  in  their  analysis  with  regard  to  the  relative  rates 
of  change  between  wavenumber  and  wave  amplitude.)  A formalism 
is  presented  in  the  present  paper  that  generalizes  Whitham's 
theory  to  include  high-order  dispersive  and  diffusive  effects 
and  which  extends  Landahl’s  (1972)  modification  of  the  basic 
action  equation  for  linear  nonconservative  effects  to  fully 
nonlinear  ones. 

Two  additional  references  can  be  cited  in  the  context  of 
the  present  discussion.  Bretherton  (I969)  was  able  to  show 
how,  for  linear,  nondissipative  problems,  the  phase-averaged 
Lagrangian  density  integrated  over  modal  cross-space,  satisfies 
.Whitham’s  action  law,  thus  extending  its  applicability  to 
wave  motions  in  waveguides.  Under  the  stated  restrictions 
Brethejrton* 3 theory  was  extendable  to  higher  order,  and  in 
principle,  the  high-order  dispersive  corrections  as  previously 
discussed  are  reproducible,  Bretherton* s low-order  modal 
results  were  generalized  to  include  nonlinear  effects  by 
Hayes  (1970b),  who  in  addition  postulated  a local  conservation 

t 

law  in  the  cross-space  divergence  operator.  However,  in 
addition  to  its  restriction  to  purely  conservative  problems, 
Hayes*  formalism  relied  on  certain  "well-known"  variational 
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identities  which  cannot  be  extended  to  avstems  of  arbitrary 
order.  In  the  analysis  to  be  presented,  both  Bretherton's 
and  Hayes'  generalizations  of  Whitham's  ideas  too  are 
extended  so  that  they  appear  as  specific  subsets  of  a broader 
unifying  theory.  The  mathematical  approach  taken  considers 
a general  variational  principle  of  arbitrary  order  in  which  the 
explicit  presence  of  modal,  propagation  space  and  time 
coordinates  is  accounted  for,  in  addition  to  that  of  the  work 
function  describing  general  nonconservative  effects.  Slowly 
varying  wave  solutions  are  sought  in  the  form  of  a Luke  ( ) 
transformation  modified  for  modal  dependence,  and  generalized 
phase  and  amplitude  equations  are  derived  from  existence  and 
secularity  arguments,  A number  of  specific  applications  are 
made,  leading  to  a number  of  formulas  useful  in  direct 
applications  of  the  general  theory,  and  in  the  case  of  weakly 
nonlinear  waves,  the  structure  of  the  high-order  dispersive 
and  diffusive  terms  is  determined  from  a pseudo -analytic 
continuation  method.  This  general  theory  describing  the 
dynamics  of  kinematic  waves  completes  the  first  part  of  the 
paper,  and  in  the  second  part,  some  new  ideas  in  hydrodynamic 
(and  general  wave)  stability  are  discussed.  The  physical 
bases  behind  some  of  these  ideas  are  discussed  in  detail  using 
various  analogies  from  gasdynamics,  but  because  of  obvious 
difficul  tties,  detailed  calculations  are  not  pursued  in  the 
present  paper.  Some  of  the  ideas,  admittedly,  are  speculative, 
but  it  is  the  author's  hope  that  the  discussions  presented 
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herein  stimulate  further  investigation  and  study.  In  a sense, 
the  second  part  of  the  present  paper  deals  with  novel  ideas 
that  hopefully  would  be  received  with  the  same  enthusiasm 
as  were  Whitham's  (I967)  ideas  on  "elliptic  instability"  and 
wave  group  "splitting". 

The  analysis  that  follows  essentially  summarizes  the 
main  results  of  a doctoral  thesis  by  the  present  author 
(Chin,  1976)  presented  at  the  Massachusetts  Institute  of 
Technology,  Further  papers  are  planned  and  the  reader  is 
urged  to  refer  to  that  work  for  details  of  ideas  and 
calculations  briefly  discussed  here.  Among  the  subjects  planned 
for  future  discussion  are  theories  detailing  the  wave 
back- interact ion  on  inviscid  shear  flows,  the  structure  of 
foci  that  result  from  wave  trapping  due  to  inhomogeneities,  and 
other  aspects  of  wave  stability. 
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II.  The  General  Theory.  Low-Order  Results. 

Let  U3  suppose  that  the  physical  system  under ' consideration 
is  derivable  from  the  variational  principle  as  give.^  in  Eq,.  (1), 

V,  * 

with  Lagrangian  density  L and  dissipation  functional  P.  We 
consider  variational  principles  (as  opposed  to  specific 
equations)  so  that  results  common  to  all  systems  can  be  obtained. 
Also  let  the  small  parameter  E characterize  an  appropriate  slow 
variation  which,  in  addition,  is  proportional  to  the  weak 
influence  of  nonconservative  effects. 


L + L' 


L I-  U,  ) 

L - L*  « L 


e F(u,...) 


(1) 


'MX, 


U 


(-0 


In  Eq.  (1)  sums  over  various  e<,  ^’s  are  understood,  x and 

y being  the  propagation  and  modal  coordinates,  and  t being  the 
time.  The  symbols  o<,  |g,  ^ indicate  the  respective  order  of 
partial  differentiations  over  x,  t and  y.  Further,  it  is 
assumed  that  for  F » 0 the  system  admits  propagating  waves.  We 
now  attempt  an  asymptotic  solution  in  the  form 
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U(x,i),-t)  - LTCe.X,!,  ys  t)  (2) 

X*M  T'ti  0(x,-l:)  - e''fli(X,T) 
a)(XJ)— e^-0^--2)(X.T)  K(X,T}-9,=  ix  ' 


In  the  above  I X and  T are  stretched  space-time  coordinates 
chosen  to  account  for  0(1)  variations  in  U over  large  scales, 

6 is  the  phase  function  for  periodic  behaviour,  and  K and  CU 
are  wavenumber  and  frequency  satisfying  the  condition  of  wave 
conservation.  We  shall  neglect  any  exponentially  small 
reflections  and,  in  addition,  for  brevity  neglect  high-order 
scaling  variables  and  multiple-phases  (which,  as  will  be  seen, 
can  be  easily  incorporated  into  the  present  approach).  The  key 
to  a general  analysis  is  the  definition  of  the  operator  H in 
the  space  9,X,T 


Uj,K  ’k 


(3) 


Thus,  neglecting  terras  of  0(t*),  the  Euler  equation  specified 
in  Eq,  (1)  can  be  rewritten  as 
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F(u,...) 


For  reasons  which  will  later  be  apparent,  we  express  Eq.  (4) 
in  the  alternative  form 

^ H^L)  - (U)C-|.  m\  * I = eUeF(U)  (5) 

where  L*  has  been  eliminated  using  the  identity 

- L'U,  (6) 

This  form  of  the  Euler  equation  leads  to  a convenient  derivation 
for  a "generalized  action  law"  which,  as  will  be  seen,  follows 
as  a necessary  condition  for  the  existence  of  assumed  asymptotic 
- solutions  in  the  form  given  in  Eq,  (2),  Now,  twice  application 
of  the  differential  identity 

to  Eq,  (5),  first  with  primes  denoting  d derivatives  (where 

I ^ 

g « U-  and  f = ) and  second,  with  primes  denoting 

" jy 

y-derivatives  leads  directly  to 


•0 

A 
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eU,F(U) 


(8) 


2>  L \J^*^^)e,rij  ] 

|+fl<K  L^rjj, (<+!)», X ^ , 


.•‘.j-K.  r 


•CH  f ^ ^ ^ *) 

+ eocK  2?  l^,(rH)jjH^4Wtf“Lx,(w)||^*‘V)».jM 

® { l-u-o^^«it#)9 ' ^(^2)^)^t♦#)9,|j  ■*’-}  ^ 

-o^CkJ) 


where 


i K 22  {L,g^^4l^j9Ug  +-} 


(9) 


■*■  e«c  K 22^  { L iy,c<4|.2)9,X^e  tijiCt+y-J)*.!  ) 
■*‘«^K  22  {Lyy^t^.2j9^7U5  "^3,(m+h)»,t^>9'^“1 


•4, 
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L 


We  now  expand  the  variables  L»  Lt  U and  } in  series  like 

(1/ » i E"v(/‘"*Ce.X.T)  (lo) 

and  introduce  the  resulting  expressions  in  Eq,  (8),  By  equating 
coefficients  in  the  varying  powers  of  S-  » the  leading  terms  give 


Averaging  over  phase*  we  obtain  the  "cross-space  conservation  law* 


a generalisation  of  Hayes'  (1970)  result.  The  bracketed  quantity 
above  can  be  interpreted  as  an  energy-like  density  that  varies 
slowly  in  time  t and  space  x.  Since  L and  U both  depend  on  x 
and  t,  Eq,  (12)  can  be  thought  of  as  a kind  of  modal  energy 
redistribution  law,  the  exact  nature  of  which  will  depend  on 
the  dynamics  of  the  propagating  wave.  Expansion  of  Eq,  (8)  and 
retsntioniof  0(  £ ) terms  gives 
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»4(lV/)  - l/rFCU™} 


(13) 


+ K/'CC.,  I 

ir^'^  7r‘*^  -+  I 

1 " '-T,(jt-j)||  -J 

+ ^-«  flfrt*)  ,,c«)  re)  ift«J  ^ I 

oCK  1 ~i 

')\j  ■ 

W \f  -t'*'  1 

•i « f iT^'^  +T*'^  'i  1 

-K  \ ^ ^ni-i)i}^cc(-*0*O9y  | 

i-ft'"  u"'  +1™  U‘”  V 


In  th*  above  the  bracketed  quantities  are  arranged  so  that  the 
first  contains  time  derivatives  only  and  the  second  contains 
propagation  space  derivatives.  By  introducing  an  "average 
Lagrangian"  with  the  definition 
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B(l«  (13)  takes  on  a particularly  simple  form  on  phase  and 
cross-space  averaging.  Assuming  zero  work  interaction  at  the 
boundary  (yj^»y2)#  obtain 

hC  - = M { vrF(l/’V  (15) 

where  we  have  denoted 


^ — 


(16) 


It  follows  that  Bqs.  (11)  - (15)  are  the  basic  governing  laws 

when  solutions  of  Eq.  (1)  are  sought  in  the  form  expressed 

by  Eq.  (2).  These  results  are  easily  generalized  to  higher 

dimensions  by  taking  ( ) "*■  V-*  $ ( )_  • »nd 

y y A jc 

introducing  the  wave  irrotationality  requirement  that 
Vjf  X IT  = 0 in  addition  to  "*■  Eq.  (15)  is 

precisely  Whitham's  action  law  when  F * 0 and  y-dependences  are 
suppressed.  For  such  "local"  waves*  Whitham  (1970)  has  shown 
(using  L « L(u,u^,u^))  that  wave  action  conservation  is 
"absolute"*  i.e.*  correct  to  all  orders  in  e . On  this  basis* 
Eq.  (15)*would  be  true  with  superscripts  removed.  Quite  often 
this  "absoluteness"  has  been  taken  for  granted  for  systems  of 
arbitral^  order.  But  this  move  is  not  altogether  clear*  and 


for  this  reason*  "local”  waves  are  reexamined  here. 

If  in  Bq.  (1)  y-dependencea  are  suppressed  by  setting 
y = 0,  a tedious  but  straightforward  extension  of  the  above 
analysis  shows  that*  in  general* 


(17) 


Thus*  ¥rave  action  conservation  is  not  "absolute"  for  systems  of 
arbitrary  order.  However*  when  the  "order"  of  the  system  is 
two*  it  is*  because  the  derivatives  • ®'tc,  * vanish 

identically.  In  fact  if  we  denote  by  A the  wave  action  in 
Eq.  (17)*  one  can  demonstr/ite  the  existence  of  high-order 
dispersive  corrections  like  that  modify  the  low-order  action 

equation  of  Whitham's  theory*  and  for  dissipative  systems* 
high-order  Burgers’  type  diffusive  terms  that  involve  (this 
will  be  shown  later  in  Section  III).  The  above  equation  can 
also  be  obtained  in  a more  powerful  way  by  extending  some  of 
Whitham's  ideas.  Let  us  consider  the  three -variable  variational 
principle 

XX  XX 

sff]  L(X,T,U.H(U))d0<JXclT  - efll  F(\J)  SUdeJXJT  - 0 (18) 


for  the  three-variable  function  U(  d*X,T).  Here  U and  its 
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variations  are  assumed  to  be  periodic  in  0$  and  variations  in 
0 vanish  on  the  boundary  R.  Variation  of  U gives,  after  some 
manipulation,  the  two-timed  Euler  equation,  Eq.  (^)  with  'H  = 0, 
Thus,  the  two-timed  form  of  the  Euler  equation,  Eq,  (1),  is 
precisely  the  Euler  equation  obtained  from  the  two-timed 
variational  principle.  It  follows  that  Eq,  (18)  is  exact  because 
it  contains  the  whole  expansion  and,  in  addition,  it  is  already 
the  averaged  variational  principle.  We  now  introduce  an  amplitude 
measure  "a"  so  that  the  solution  for  U(  d,X,T|£  ) depends 
explicitly  on  the  two  parameters  "a"  and  6 , Then,  since 

slfEJXJT  ' ii^{kjF(v)W8*i;jii)c(e}dXJ7  (w) 

R T?  0 

mist  hold,  **0"  variations  lead  to 

and  6 variations  lead  to  Eq.  (17).  Equations  (17)- (20)  are 
then  the  required  "generalized  action  and  phase  relations*. 

For  low-order  solutions  the  choice  U(  5»X,T,£  ) = U^*^(a,  9 ) 
where  U^*^^a,  8 ) is  the  uniform  wave  solution  to  Eq.  (1)  suffices, 
as  shown  by  Whitham  (1970),  But  for  proper  account  of  high-order 
effects,  it  is  in  general  necessary  to  assume 
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iroATiO  - 1 (21) 

5*0 

in  view  of  the  work  by  Yuen  and  Lake  (1975).  This  point  will  be 
discussed  more  fully  in  Section  III.  The  proper  choices  for 
the  functions  are,  in  general,  unclear.  This  difficulty 
is  in  addition  compounded  by  the  fact  that  variational  principles 
are  not  unique  - a concern  that  does  not  arise  in  the  low-order 
problem  because  the  phase  and  action  equations  are  unaffected. 

If  Eq.  (21)  is  chosen  correctly,  though,  the  analysis  carried 
to  the  appropriate  level  of  approximation  does  produce  high-order 
dispersive  terms  (as  are  modelled  in  recent  studies  by  a certain 
nonlineEU*  Schroedinger  equation)  and,  in  addition,  high-order 
diffusion  terms  that  can  be  important  near  kinematic  shocks. 

Now  to  loading  order,  the  evolutionary  equatio7«»  are 

. If-  (22) 

L‘"  - o u(e,x,Tit)  “ U‘"(a,0) 

and  it  is  known  (Whltham,  19^7)  that  the  physical  properties  of 
the  solution  vary  markedly  according  to  the  type  of  the  governing 
system.  For  example,  Benjamin- Feir  or  Whitham-type  modulational 
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Instabilities  are  encountered  for  elliptic  systems.  For 
hyperbolic  systems,  the  perturbation  "splits"  and  propagates 
at  two  different  speeds.  And  in  the  linear  (parabolic)  case 
the  modulation  equations  uncouple,  and  the  wavenumber  field 
is  determined  independently  of  the  wave  amplitude.  In  this 
case,  Eq,  (22)  shows  that  the  concept  of  group  velocity  ^ still 
relevant  to  weakly  nonconservative  processes.  It  is  still 
defined  by 

C " ” 5ai/aK  (23) 

since  the  phase  relation  can  be  written  in  the  form 

a a^jiS  (Cu,K,X,T)  a 0,  To  this  order,  the  kinematics  are 
independent  of  the  functional  P,  agreeing  with  the  fact  weakly 
nonconservative  effects  affect  the  phase  only  to  0(£^).  (A  more 
■ explicit  formula  showing  the  general  nature  of  high-order 
dispersive  and  diffusive  corrections  to  Eq,  (22)  will  be  given 
later  in  Section  III,)  Wave  action  is  still  defined,  and  again, 
in  terms  of  the  associated  conservative  system, 

A - t:  (24) 

but  it  is  dissipated,  in  general,  with  strength  UgF  • The  rays 
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»re  still  determined  by  the  "real  group  velocity",  where  we  have 
avoided  the  inconsistencies  of  the  complex- frequency  approach. 
These  restilts  apply  to  fully  nonlinear  waves  as  well. 

In  practice  it  is  sufficient  to  choose  for  the  uniform 
wave  solution  to  the  associated  conservative  system.  However 
the  presence  of  the  right-side  integral  of  Eq,  (22)  in  nonlinear 
systems  may  possibly  alter  the  structure  of  the  fundamental 
conservative  solution  in  a significant  way.  This  point  is 
especially  important  in  bifurcational  and  transitional  studies. 
We  may  also  note  that  by  considering  dl/dT  and  dl/dX  instead 
of  "bl/bB  in  Eq,  (6),  one  obtains  for  local  or  modal  waves  the 
following  leading  order  results, 

#r(wL:-r)  - a^Cwi;:')  » -t;'  (25) 

- h ( kC)  - < l'k"  -r)  - L'x  - K (25) 

These  identify  Uq  as  the  appropriate  weighting  function  for  F, 
The  above  equations  are  equivalent  to  Eq.  (22),  except  that 
different  energy-liJce  norms  are  now  examined  along  different 
rays.  As  for  Eq.  (22),  the  superscript  zeroes  may  be  removed 
, for  second  order  systems,  but  for  high-order  systems,  additional 

j terms  must  be  incorporated.  If  we  identify  Eq.  (1)  with 

Hamilton's  principle,  we  see  that  the  above  represent  energy 
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and  momentum  laws  respectively.  This  follows  for,  in  a 
mathematical  sense,  energy  is  that  quantity  conserved  when  the 
corresponding  variational  principle  is  invariant  with  respect 
to  translations  in  time.  The  conservation  laws  for  momentum 
and  action,  for  example,  follow  from  considering  translations 
in  space  and  in  phase.  The  right  sides  of  Eqs,  (25)  and  (26) 
show  explicitly  the  role  of  inhomogeniety  and  dissipation. 

For  exajnple,  the  inhomogeniety  can  be  prescribed  (as  would  be 
the  case  in  linear  theory)  or  can  be  the  result  of  wave- induced 
modifications  to  a given  mean  state.  In  linear  tlieory  the 
velocities  for  action,  momentum,  energy  and  wavenumber  are 
identical,  of  course,  since  the  average  Lagrangian  is  homogeneous 
in  a^. 

It  is  important  to  recapitulate  several  ideas  fundamental 
to  the  above  analysis.  The  first  is  that  no  assumptions  have 
been  made  regarding  the  orders  of  the  conservative  and 
nonconservative  operators  in  Eq,  (1),  the  variability  of  their 
coefficients,  or  their  linearity  or  nonlinearity.  Secondly, 
the  general  expansions  introduced  in  Eq,  (2)  include  a strong 
modal  dependence,  as  opposed  to  those  used  by  Luke  ( ) or 
Whitham  (1970),  One  observes  that  when  Eq,  (10)  is  introduced 
in  Eq,  (8)  with  y=  0,  the  left  side  of  the  resulting  0( £ ) 
equation  contains  functions  of  both  U^^^and  while  the  right 
side  contains  exclusively,  A solution  uniformly  valid  in 
0 requires  U,  and  hence  each  to  be  periodic  in  0 with 
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period  ZX  , This  uniform  validity  is  possible  only  if  the 
integral  of  the  right  side  with  respect  to  6 over  a period 
is  zero.  This,  together  with  the  definition  given  in  Eq,  (14), 
leads  to  the  action  law  (Eq.  (15))  directly  as  straightforward 
manipulations  show.  In  other  words,  wave  action  conservation 
follows  as  a secular  condition  (in  the  sense  of  Cole  (I968)) 
when  solutions  of  Eq,  (1)  are  sought  in  the  form  specified  by 
Eq,  (2),  For  the  treatment  of  modal  effects  eind  high-order 
corrections,  secularity  alone  is  not  sufficient,  but  phase 
averages  can  still  be  conveniently  carried  out  via  the 
two-timing  introduced  in  Eq,  (2), 

It  is  crucial  that  Eq,  (22)  holds  for  equations  of  any 
order,  but  that  the  order  of  the  system  becomes  important  only 
when  higher  order  effects  are  to  be  calculated.  For 
second-order  systems,  however,  Eq,  (22)  is  exact  because 
derivatives  of  frequency  and  wavenumber  vanish  identically,  and 
superscripts  may  be  dropped.  Thus,  Whitham's  (1970)  comments 
on  the  absoluteness  of  action  conservation  (in  nondissipative 
media)  apply  only  to  these  degenerate  cases.  Action  conservation 
here  turns  out  to  be  the  appropriate  adiabatic  invariant  for 
linear  or  nonlinear  wave  propagation  in  inhomogeneous  dissipative 
systems.  In  nonlinear  processes,  this  is  especially  important 
because  energy  cannot  play  too  fundamental  a role  - it  is  easily 
transferred  between  component  frequencies.  For  solutions  to 
practical  problems  it  is  necessary  to  amend  Whitham's  (1970) 
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algorithm.  In  linear  problems,  for  example,  the  uniform  wave 
solution  U^*^=  as  in  8 led  to  an  average  Lagrangian  hmogeneous 
in  a^  I this  is  true  only  in  the  lowest  approximatio  v.  Since 
•a"  is  really  a function  of  the  slow  variables  X and  T,  0(£^) 

corrections  to  the  average  Lagrangian  must  be  incorporated 

/ 

that  will  involve  space-time  derivatives  of  wavenumber  and 
amplitude,  0( £ ) effects  having  vanished  by  orthogonality.  The 
existence  of  such  high-order  terms  suggests  a new  treatment  for 
"kinematic  shocks".  There  thus  exist  high-order  diffusion 
terms  that  can  possibly  smooth  out  discontinuities  that  appear 
in  considering  the  low-order  solution  alone  or,  for  example, 
high-order  dispersive  terms  that  may  suppress  shock  formation. 

The  dispersive  corrections  mentioned  here  (as  will  be  pointed 
out  more  clearly  later)  resolve  the  inconsistencies  pointed  out 
by  Chu  and  Mei  (1970)  regarding  the  apparent  loss  of  certain 
high-order  modulation  terms  in  the  Lagrangian  formulation  by 
showing  that  such  terms  are  accounted  for  by  the  variational 
method.  It  is  only  that  through  the  above  inconsistency  that 
the  discrepancy  arises. 

Our  approach  also  bypasses  the  conventional  use  of  complex 
frequencies  in  treating  dissipative  problems.  For  example,  it 
is  not  necessary  to  determine  the  group  velocity  from  the  real 
part  of  a,  complex  frequency,  and  then,  to  introduce  heuristically 
into  the  energy  equation  a damping  based  on  the  imaginary 
frequency.  In  many  nonlinear  problems  where  the  use  of  complex 
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frequencies  is  inappropriate,  Eq.  (22)  seems  to  provide  the 
only  practical  alternative.  Both  local  and  modal  w-ves,  of 

course,  were  considered  under  the  same  formalism,  t lere  is 

\ * 

still  another  interesting  point  to  be  noted.  For  nonlinear 
waves  we  can  eliminate  amplitude  from  the  argument  space  of 

i 

and  using  the  dispersion  relation.  The  action  equation 
then  becomes  an  eikonal  one,  from  which  an  equation  for  a 
perturbation  phase  could  be  developed.  For  this  system  it  turns 
out  that  boundary  conditions  cannot  be  satisfied  if  spatial 
growth  is  not  assumed,  for  initial  values  applied  at  x = 0, 

On  the  other  hand,  initial  values  applied  at  t = 0 require 
temporal  growth.  This  extends  in  a simple  way  Gaster's  (1965) 
result  that  Orr-Somraerfeld  waves  (generated  from  a vibrating 
ribbon  fixed  in  space)  on  a boundary  layer  grow  spatially.  It 

also  reduces  to  Taylor's  (1962)  results  for  temporally  growing 

I 

waves. 

In  applying  Eq,  (15)  to  specific  problems,  the  choice  of 
the  basic  solution  (or  Ansatz)  is  crucial.  From  the  linear 
results  to  be  presented  (Example  1),  the  uniform  solution  of 
the  associated  conservative  system  appears  to  suffice  as  a 
general  rule  because  the  dissipation  is  weak.  This  model  can 
be  extended  to  nonlinear  and  modal  problems  as  well.  For 
example,  ,if  Luke's  (1966)  nonlinear  Klein-Gordon  equation  were 
modified  to  include  dissipation,  one  would  evaluate  Eq,  (15) 
by  the  uniform  conservative  solution  he  obtained.  Finally 
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not*  how  the  conserved  density  In  Eq.  (12)  takes  on  a much 
more  complex  structure  than  the  analysis  of  Hayes  (1970)  would 
suggest,  but  reduces  to  his  in  the  appropriate  limit.  This 
comment  applies  to  the  modal  law  Eq,  (13)  es  well,  noting  that 
neither  Bretherton* s (I969)  nor  Hayes'  theories  are  valid  for 
dissipative  wave  flows.  The  evolutionary  equations  derived 
have  a two-fold  use.  First,  they  can  be  used  to  describe 
the  mathematical  properties  of  wave-like  solutions  to  specific 
equations,  and  this  is  the  usual  application.  The  second  is  to 
make  certain  physical  assumptions  about  the  functional  form  of 
the  Lagrangian,  for  example,  that  it  represents  Hamilton's 
principle,  and  then,  obtain  alternative  descriptions  for  physical 
conservation  laws.  This  approach  is  pursued  later  in  an  example. 
One  of  the  aims  of  the  present  paper  is  a systematic  re-derivation 
and  extension  of  both  existing  high  and  low-order  theories  of 
wave  motion  from  one  general,  self-consistent  approach.  Before 
considering  in  detail  the  high-order  theory,  we  discuss  several 
applications  of  the  low-order  theory. 
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Example  1«  Linear  "Local"  Waves. 

First,  one  can  obtain  specialized  results  for  rather 
general  linear  systems.  Consider  "local  waves",  that  is,  waves 
without  a modal  cross-structure.  We  examine 
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where  the  sums  are  taken  over  combinations  of  + ^ = even 
and  r+  y = odd.  The  coefficients  a_^  , and  b ^ are  assumed  to 

*‘tp  ’ #• 

be  weakly  dependent  on  x,t  and  to  be  sufficiently  "well-behaved". 
The  conventional  approach  assumes  a solution  proportional  to 
exp  i(Kx-Wt),  where  is  complex  and  uJi/w^  ^ Qi  t ). 

Neglection  of  0(£^)  terms  in  the  complex  dispersion  relation  gives 
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Equating  real  and  imaginary  parts  we  find  the  relations 
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J-Z  l’’*”  <^r'‘ 


0(1) 


(29) 


(30) 
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Solutions  corresponding  to  the  neglected  roots  of  Eq,  (28)  have 
been  discarded  because  they  are  associated  with  rapidly  varying, 
or  alt t. -natively,  highly  damped  solutions.  Thus,  if  we  are 
looking  for  slowly  varying  wavetrains  in  which  dissipation  plays 
a modest  role,  Eq,  (28)  must  suffice,  Eq,  (29)  can  be  interpreted 

t 

as  a dispersion  relation,  consistent  with  that  in  Eq,  (22)  in 
the  sense  that  it  is  not  affected  by  weakly  nonconservative 
effects.  On  the  other  hand,  Eq,  (30)  provides  a compact 
expression  for  the  "exponential  growth  rate".  Now,  assume  that 
the  left  side  of  Eq,  (2?)  is  derivable  from  a Lagrangian  density, 
so  that 


- L 


u.'^  9 X 


iL 


(31) 


XX 


i 
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since  the  term  aQ  q(X,T)  does  not  alter  the  analysis  we  make, 
without  loss  of  generalization,  the  normalization  a^  ^ = 1, 

Then  the  Lagrangian  must  take  the  form  L = iu^  , this 

"fixed"  term  enabling  us  to  keep  track  of  the  entire  conservative 
operator  throughout  the  averaging.  From  this  it  follows  that 


(32) 


when  L is  evaluated  with  the  uniform  "conservative"  wave 
■ solution  a(X,T)8in8  . Thus  unlike  conventional  methods 
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in  which  diasipation  is  separately  accounted  for  in  a complex 
phase,  the  "a"  in  our  approach  contains  both  conservative  and 
slightly  nonconservative  effects.  Evaluating  the  right  side 
of  Eq,  (27)  with  this  solution  gives,  successively, 


■2L‘;’ 


-ri.)  ?5:b,_r*‘"KV 


(33) 
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T 


In  terms  of  the  wave  action  A = and  the  group  velocity 
C = we  have 

- 3^  ^ r>  ‘ 

Which  describes  the  long  term  effects  of  locally  weak  dissipation. 
It  is  essential  here  that  the  basic  group  velocity  (which 
determines  the  kinematics)  remains,  to  lowest  order,  unaffected 
by  dissipation.  The  above  equation  has  been  used  by  Landahl  (1972) 
and  by  Davey  (1972),  Complex  frequencies,  again,  are  unnecessary 

t 

in  the  present  approachi  the  purpose  of  Eq,  (3^)  io  to  relate  A 
to  the  more  common  "U)j", 
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Example  2,  Linear  "Modal"  Waves, 

Modal  waves  are  those  that  involve  a cross-structure 
orthogonal  to  the  propagation  direction.  A simple  model  is 


* I 


^ U 


(35) 


where  the  left  side  is  conservative  and  the  right  side  is 
dissipative.  Now  define  the  auxiliary  system 


(36) 


having  the  eigensolution  ^ ®*P  i(Kx-u»t)  with  mode 

number  "s",  where  9^3(7)  is  determined  from  the  solution  to 


and  its  appropriate  homogeneous  boundary  conditions.  The  right 
hand  sides  of  the  above  may  be  dropped  if  the  resulting  undamped 
system  is  free  of  critical  layers.  Then,  the  extension  to 
Whitham* Si  basic  algorithm  suggests  an  Ansata  of  the  form 


U;’'(S,X.T,J)  - «s(X,T;9’(3)im9 


(38) 


This  type  of  assumption  is  reminiscent  of  Stuart's  (1958) 
"shape  assumption",  and  appears  to  be  justifiable  of  physical 
grounds  (provided  that  the  medium  is,  indeed,  slowly .varying) . 
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Example  3,  Waves  in  Moving  Med  la t Irrotatlonal  Mean  Flowa, 

The  evolutionary  equations  as  given  in  Eq,  (22/  do  not 
impose  any  restrictions  on  the  background  medium,  except  that 
it  is  slowly  varying.  We  will  consider  flu id -dynamical 
problems  in  which  waves  ride  on  slowly  varying  mean  currents 
that  vary  slowly  in  time  and  in  the  propagation  direction  only. 
Both  wave  and  mean  flow  interact  nonlinearly  through  a radiation 
stress  that  transfers  energy  from  mean  flow  to  wave  (and  vice 
versa).  In  addition,  the  slow  variation  of  the  medium  introduces 
small  local  changes  to  the  dispersive  properties  characteristic 
of  the  wave.  The  cumulative  effects  of  these  changes  may  be 
significant. 

Let  us  assume  that  the  waves  under  consideration  propagate 

in  a direction  parallel  to  the  mean  flow  U ~ U(X,T),  This  base 

\ 

flow  is  irrotational  because  there  is  no  dependence  on  the 
shear  coordinate  y.  Then  if  superscript  primes  refer  to  a 
coordinate  system  S'  moving  with  the  mean  flow  U and  unprimed 
symbols  refer  to  ground-fixed  coordinates  S,  the  principle  of 
Galilean  invariance  implies,  to  leading  order, 

U)-UK,  K)  (39) 

I 

What  enters  in  Eq,  (39)  is  the  expected  Doppler  shift,  as  can 
be  derived  by  applying  Eq,  (2)  and  its  primed  counterpart  to 
any  system  of  equations  written  for  both  S and  S'  (Chin,  1976), 
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The  key  point  in  the  analysis  that  follows  is  a re- interpretation 

of  the  fictitious  forces  experienced  in  S’  as  equi\ ^lent  Reynolds 

stress  effects  as  observed  from  S,  This  principle  .vaa  first 

enunciated  by  Whitham  (1962)  in  the  context  of  gravity  waves. 

It  may  have  been  suggested  by  the  fact  that  the  flu id -dynamical 

equations  linearized  about  a constant  U are  identical  with  the 

same  equations  written  in  moving  coordinates,  but  the  analogy 

is  not  quite  complete.  With  the  above  interpretation  simple 

results  are  available  on  defining  energy  densities  E' 

and  E (which  are  valid  for  fully  nonlinear  waves),  as 

10 

suggested  by  Eq,  (25),  Now  because  wave  action  is  invariant 
under  Galilean  transformation,  that  is,  -C  , we  can 
write  (to  leading  order)  (E’ )/ OU'  ~ (E+£)/cu  - • 

For  linear  systems,  Ji  = X,'  - ® since  the  difference  between 
phase-averaged  kinetic  and  potential  energy  densities  vanishes. 

If  the  resulting  expression  for  is  substituted  in  Eq,  (22) 
and  Eq,  (33)  is  used  for  the  dissipation  term,  we  are  led  to 
(on  dropping  primes) 


;)E=  Ej^hl 


un.  - Kn.u, 


+ 


where  wet  have  denoted  jQ(K,X,T)  = u/*  and  set  CO  s U(X,T)K'*‘/I 
in  the  wave  conservation  law.  The  group  velocity  here  is 
simply  Cg  = It  is  easy  to  see  that  the  right-side 

I ' term  in  Eq,  (40),  aside  from  the  Wi  , represents  the  effects 

i 
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of  radiation  stress.  For  example,  progressive  waves  on  deep 
water  satisfy  Si  - gK  + T^^K^/^  where  g is  the  accele 'ation  due 
to  gravity,  ?^  is  the  surface  tension,  and  ^ is  the  .fluid 
density.  Since  = Sl^  = 0 it  follows  that  the  radiation 
stress  R = (1+ 3Z)E/(2 -*■  2Z)  where  Z = For  pure  gravity 

waves  (Z  = 0)  R=iE,  and  for  pure  capillary  waves  (Z  = oo  ) 

R = 3E/2,  reproducing  the  results  of  Longuet-Higgins  and 
Stewart  (1964)  in  a simple  way.  For  gravity  waves  in  water  of 
finite  depth  h(X,T),  the  intrinsic  frequency  is  SI  - (gKtanh  Kh)^. 
In  this  case  Sl^^  and  SL^  are  nonzero.  Expanding  these  terms 
out  and  simplifying  with  h^  + (hU).^  = 0 produces  the  interaction 
term  (J-+ 2Kh/sinh  2Kh)E  in  Eq.  (40),  which  again  agrees  with  the 
results  of  detailed  perturbation  methods.  These  results  are 
not  altogether  surprising  in  view  of  the  physical  ideas  underlying 
Eq,  (39),  Eq,  (40),  of  course,  is  valid  for  non-fluid-dynamical 
problems  as  well.  We  simply  set  U = 0 to  arrive  at  the  interaction 
term  E SL^/St  . 

In  contrast  spatial,  as  opposed  to  temporal,  inhomogenieties 
are  responsible  for  changes  in  wave  momentum.  The  mathematical 
basis  for  this  is  implied  by  Noether's  theorem  (Courant  and 
Hilbert, ■ 1953).  A wave  momentum  law  analogous  to  Eq,  (40) 
can  be  easily  derived.  For  linear  systems,  this  is  accomplished 

I 

by  noting  that  E = /lA  and  M = KA  for  the  wave  quantities. 

Thus  E/y2  * M/K  which,  together  with  action  and  wave  conservation, 
leads  to 
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(41) 


It  is  also  possible  to  extend  Eq.  (40)  to  completely  nonlinear 
flows.  This  is  accomplished  by  applying  Eq.  (39)  to  Eq,  (25)i 

I 

with  the  result  that 


(42) 


where  now  E is  the  wave  energy  density.  One  may  have 

anticipated  Eq.  (42)  in  retrospect.  For  example,  the  right-side 
terms  are  for  radiation  stress,  since  -£>  is  a morae.itum 

flux  (see  Eq,  (26)).  The  momentum  analogue  to  Eq,  (42)  is 
just  as  easily  derived  by  combining  Eqs,  (39)  and  (26),  The 
result  is  simply 


H 


i^)M  - 4 -M 


dx 


(43) 


The  above  laws  are,  of  course,  coupled  to  the  dynamics  of  the 
mean  flow.  For  further  details,  the  reader  is  referred  to 
Whitham  (I967),  Three  propagation  velocities,  those  for  action, 

I 

energy  and  momentum,  can  be  defined.  They  are,  respectively. 
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(1*4) 

c“’  ’ V - 

C'”'  ‘U*U-Ki,yu^  '■:■ 

In  the  linear  case  these  reduce  to  U + fl^  since  the  Lagrangian 

/V 

vanishes.  Note  also  that  the  nonlinear  wave  kinetic  energy  K 
satisfies  exactly  K = ) = ^ co  A = i(i^/K)(KA)  = iCpM  , 

where  is  the  phase  velocity.  This  identity  is  based  on  E 
being  the  sum  of  kinetic  and  potential  energies  and  ^ being  the 
difference.  It  was  first  given  by  Levi-Civita  (1924)  for  finite 
amplitude  gravity  waves  in  water  of  finite  depth.  But  its 
applicability  seems  to  extend  to  much  more  general  systems.  In 
the  linear  case,  energy  equipart ition  implies  E = 2 K and  so, 

E = CpM,  a result  familiar  in  surface  waves.  ^ 

General  results  like  Eqs.  (42)  and  (43)  were  made  possible 
by  the  assumed  dependence  of  U on  X,T  only,  that  is,  there  is  no 
y-dependence.  Thus  the  dispersive  properties  of  the  wave  are 
unmodified  by  the  presence  of  the  mean  flow  except  to  the  extent 
of  the  expected  Doppler-shift , For  many  problems,  however,  the 
assumption  U « U(X,T)  is  not  very  realistic.  In  deep  water,  for 
example,  the  nonlinear  coupling  between  wave  and  mean  flow 
disappears  (Whithara,  196?)  so  that,  effectively,  the  mean  flow 

I 

is  unaffected  by  wave  growth  or  decay  as  induced  by  radiation 
stresses.  Surfaces  changes  in  U can  only  be  accounted  for  by 
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incorporating  a y-dependence.  Obviously,  the  above  approach 
using  Eq.  (39)  is  no  longer  applicable,  and  a theory  for 
the  y-coupling  must  be  separately  developed.  This*  'ill  be 
outlined  in  the  next  example. 

At  this  point  we  can  make  some  remarks  on  wave  packet 

/ 

growth.  Usually  stability  theory  centers  on  the  idea  of  temporal 
growth  or  decay,  for  example,  the  time  behaviour  of  a normal 
mode  is  studied  in  parallel  flow  models.  Wave-like  disturbances, 
however,  propagate  downstream  in  reality,  and  pass  through 
various  regions  which  may  be  stable  or  unstable  on  the  basis  of 
normal  mode  analysis.  Not  only  is  the  space-time  evolution  of 
the  wave  energy  density  important,  but  its  integral,  the  overall 
wave  packet  growth  in  defining  stability.  In  general,  the 
stability  (or  instability)  of  one  does  not  imply  the  same  for 
the  other.  Consider  the  integral  2!.(t)  of  E(x,t)  over  dx 
between  two  rays  x^(t)  and  X2(t).  is  a energy  in 

some  sense,  with  time  playing  a parametric  role.  Differentiation 
leads  to 


^E(x,t) 


cJx, 


dXt 


or, 


T 


39 


If  q is  identified  with  the  group  velocity  , then 

2!  = d2 /dt  is  simply  the  time  rate  of  change  of  total  energy 
between  two.  rays.  If  we  introduce  Eq.  (^0)  we  find  that 


d-t 


E(x,{) 

%Ci) 


JL 


(45) 


which,  in  non-fluid-dynamical  problems,  takes  on  U = 0,  If 
the  inhomogeniety  and  the  dissipation  vsmish,  Eq.  (45)  shows 
that  the  energy  between  two  rays  is  a constant  equal  to  the 
initial  energy.  In  general,  the  coefficient  of  E in  Eq,  (45) 
is  a complicated  function  of  K,  X and  T that  is  not  known 
until  the  kinematics  are  completely  determined.  However  if 
j2  t = = tt/,-  = 0 ( as  in  inviscid  deep-water  wave  problems) 

Eq.  (45)  automatically  implies  wave  energy  growth  whenever  < 0 
since  K Hif/H.  > 0. 

Consider,  for  example,  the  quantity 


a 


ji 


for  gravity  waves.  If  Q > 0 everywhere,  instability  is  implied, 
but  if  Q ^ 0 there  is  total  stability.  The  energy  of  a wave 
packet  thus  grows  or  decays  if,  relative  to  the  wave,  is 
negative  or  positive,  a fact  in  agreement  with  observation. 

Flow  inhomogenieties  play  a dual  role  in  affecting  both  energy  j 
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density  and  total  energy.  For  example,  an  energy  singularity 
can  develop  within  a wave  packet  while  the  total  energy  decreases. 
Now,  let  Qjnax  ^min  maximum  and  minimum  values  of 


. Q In  ixyX^) , so  that 


Xj(t) 

EQ  dx  ^ Q 
%Ci) 


i 


XfCtJ 


Hence,  substitution  in  Eq.  (45)  shows  that  sind  can  be 

interpreted  as  upper  and  lower  bounds  for  total  energy 
amplification.  In  linear  theory  these  bounds  are  prescribed  in 
the  sense  that  U(X,T)  is  specified  beforehand.  Since  U is  slowly 
varying  the  growth  rates  Q are  0(t  ),  however,  events  occurring 
within  the  wave  packet  can  be  quite  sudden. 

We  can  also  examine  wave  momentum.  The  total  wave  momentum 

14 

= / M(x,t)  dx  following  a wave  group  is  governed,  in 

linear  theory,  by  the  equation 


dt 


%C-tJ 


(46) 


where  x,  ,(t)  are  determined  by  integration  of  the  wave 
conservation  law.  For  simplicity  consider  waves  in  inviscid, 
deep  water,  where  ^ 0,  Then,  the  inequality 


follows,  identifying  the  maximum  and  minimum  of  as  the 
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upt>er  and  lower  bounds  for  the  momentum  amplification.  For  the 
wave  action  density  A,  one  has  d/dt  / A(x,t)  dx  = 0 in  a 
conservative  system.  For  linear  or  nonlinear  problems,  we  have 

c4  * / A(X,t)Jx  'J  A(X.0)Jx  ‘J  (U7) 

where  is  the  total  Initial  action,  with  dissipation, 

Eq,  (47)  is  no  longer  true,  but  rough  estimates  for  upper  and 
lower  bounds  on  the  decay  rates  are  easy  to  obtain.  Consider 
the  case  - -^a  “ denote 

Kj,-  Max  K(x,0)  > Kjjj=  Min  K(x,0)>0  and  Jfj^=Max( -U^^)  > )(jj=Min(-Ujj) 
it  is  clear,  from  wave  conservation,  that  along  all  rays 

Kji^e  ^ K Kj^e  ” ,‘  Since  = -ZiP  K in  deep  water,  where 
"D  is  the  kinematic  viscosity,  it  follows  that 


A(^f-t)cJx  ~ ij  UJiAh  * dx 


(48) 


and  hence,  for  any  t = t , that 


Thus,  bounds  for  the  total  action  decay  (or  growth)  can  be 
readily  formulated  in  terms  of  initial  conditions  and  streamwise 


parameters.  In  the  linear  results  given  above  it  is  possible 
for  the  total  energy  or  momentum  of  a wave  packet  to  grow 
without  bound.  In  this  case  the  effect  of  nonlinearity  must 
be  included  by,  say,  expanding  the  average  Lagrangian  in  powers 
of  wave  action.  For  further  details,  the  reader  is  referred 
to  Chin  (1976). 
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Example  Waves  in  Moving  Wedlai  Mean  Shear  Flows, 

The  previous  example  dealt  with  mean  flows  constant  across 
the  modal  cross-space.  Because  the  local  mean  speed  was  unique 
over  y we  invoked  the  principle  of  Galilean  invariance,  and 
this  led  to  general  results.  This  invariance  is,  however, 
inapplicable  to  the  case  of  mean  shear  flows.  At  first,  the 
situation  appeared  to  be  additionally  complicated  by  the  lack  of 
a suitable  variational  formulation  for  the  inviscid  Rayleigh 
equation.  This  latter  point  was  dismissed  for,  although  a 
variational  principle  does  not  exist  in  the  linearized  problem, 
it  does  for  the  complete  nonlinear  one.  With  this  settled,  we 
opted  to  equivalently  average  over  phase  and  modal  cross-space 
directly  the  fluid -dynamical  conservation  laws  for  total  mass, 
total  momentum  and  total  energy.  The  formal  procedure  parallels 


Whithara’s  (1962)  treatment  for  irrotational  gravity  waves  (of 
energy  density  E and  wavenumber  K)  in  water  of  finite  depth 
- h(X,T)  and  speed  U(X,T),  the  latter  two  variables  of  which,  on 
account  of  0(E)  wave-induced  effects,  cannot  be  assumed  as  known 
a priori,  E and  K in  addition  are  taken  as  slowly  varying 
functions  of  x,t  and  the  waves  are  assumed  to  be  locally 
sinusoidal  in  the  averaging.  The  resulting  modulation  equations 
(although  E terms  are  not  explicitly  present)  describe  the 


nonlinear  wave  back- interaction  on  the  mean  state  over  scales 

I 

large  compared  to  a typical  period  or  wavelength.  The  consequences 
of  the  theory  can  be  found  in  Whitham  (1967).  | 
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Our  theory I however,  departs  from  Whithara’a  by  recognizing 
that  the  U here,  in  addition  to  experiencing  "magnitude  changes" 
of  the  Whitham  type,  experiences  "shape  changes"  due  to  the  action 
of  Reynolds  stresses.  This  consideration  does  not  arise  in  the 
irrotational  case  by  virtue  of  Kelvin’s  theoremi  an  irrotational 
flow  stays  irrotational.  What  is  needed  to  extend  the  above 
ideas  is  thus  a suitable  (phase-independent)  Ansatz  for  the  mean 
flow  in  the  phase-modal  averaging,  so  that  within  the  limitations 
imposed  by  WKB  theory,  the  simultaneous  wave  back- interaction 
process  is  properly  described,  A simple  model  can  be  constructed 
by  observing  that  waves  truly  slowly  varying  in  space  can  be 
described  locally  by  the  usual  normal  mode  theory.  The  given 
laminar  flow  U^(y)  corrected  for  wave-induced  distortions 
of  0(E)  as  introduced  by  the  Rayleigh  equation  then  can  be  put 
in  the  form  (Chin,  1976) 

= U/u)  + E 3^^  ‘ (50) 

where  9^jj(y,K)  is  the  Rayleigh  eigenfunction  suitably 
normalized,  and  E includes  the  initially  exponential  temporal 
growth.  This  result  agrees  with  that  in  Landahl  (/973)» 

The  key  ^dea  is  to  assume  a slowly  varying  mean  flow  in  the  form 


(51) 


^5 


where  is  a reference  speed.  Thus,  in  this  model,  the  local 
normal  mode  solution  for  the  corrected  mean  flow  determines  the 
form  of  a large-scale  (wave-dependent)  shape  structure  convected 
with  the  mean  flow.  This  assumption  is  physically  sound  provided 
the  waves  are  slowly  varying.  The  above  Ansatz  reduces  to 
Whithara's  when  = constant,  since  U^"  = 0,  Second-order  changes 
in  mean  pressure,  of  course,  also  must  be  accounted  for  in  the 
averaging.  When  the  final  averaged  equations  are  combined  with 
the  kinematic  requirement  for  wave  conservation  we  arrive  at  a 
system  of  modulation  equations  which,  in  Whitham's  (I967)  sense, 
can  be  hyperbolic  or  elliptic  depending  on  the  value  of  Kh, 

These  imply  a splitting  of  the  initial  wave  packet  or  a type  of 
modulational  instability.  However,  the  new  theory  generalizes 
Whitham's  by  the  inclusion  of  velocity  profile  curvature,  and 
this  introduces  a host  of  new  effects.  The  new  theory  can  be 
considered  as  complementing  the  viscous  one  of  Stewartson  and 
Stuart  (1971)  hy  dealijig  with  essentially  inviscid  flows  that 
are  not  in  "near-equilibrium",  A detailed  treatment  for  surface 
waves  and  channel/jet  - type  flows  is  planned  in  a separate 
paper,  but  a preliminary  account  of  the  main  results  is  available 
in  Chin  (1976). 


_ 
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Example  5»  Taylor’s  Waves. 

The  Lagrangian  theory  leads  to  simple  generalizations 
for  "Taylor,  waves."  The  effect  of  nonuniform  currents  with 
constant  horizontal  divergence  on  a train  of  gravity  waves 
whose  amplitude  and  length  are  constant  in  space  but  variable 
in  time  was  first  discussed  by  G.I.  Taylor  (1962),  Such 
space-independent  waves  can  be  expected  to  exist  for  constant 
= o<  if  there  are  otherwise  no  other  inhoraogenieties  in  the 
medium,  and  Taylor  has  shown  how  these  can  be  experimentally 
realized.  We  thus  consider  "deep-water"  problems.  If  we  denote 
by  A the  wave  action,  the  nonlinear  governing  equations  on 
taking  K = A = 0 become 


“ -o<  (52) 

(53) 

which  integrate  to 

K *■  K,  e'"^  (54) 

A'A,  ef'"^  (55) 
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where  is  the  initial  wavenumber  and  is  the  initial  wave 

action.  Now  A = A^K/K^  and  we  can  expand  the  wave  energy  as 
E = Af(K) iA^g(K).  Thus  the  nonlinear  solution  for  energy 
becomes 

E “ + (56) 


X A 3 

For  deep-water  gravity  waves,  f(K)  ” Sq  K g(K)  - K 

whore  g^  is  the  acceleration  due  to  gravity,  as  can  be  shown 
from  the  results  of  Whitham  (196?).  Eq.  (56)  specializes  to 


(57) 

0 


-ja 

For  linear  waves  E = i ^ ^o*^^  implies  that  a ^ A , reproducing 
Taylor's  result.  For  capillary  waves  with  surface  tension  T, 
f = K^/2  g - -kVs  (Chin,  1976).  In  this  case. 


(58) 


-5/2 

which,  for  small  amplitudes,  gives  E ^ ^ . 

We  can  also  formulate  the  above  ideas  somewhat  differently. 
The  energy  equation  as  given  in  Eq.  (42)  is  simplified  by 
expanding  the  average  Lagrangian  in  powers  of  A.  If  wo  assume 
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that  a;  = UK  + f Ag,  the  equation  satisfied  by  (the  nonlinear 
definition  of)  E is,  to  O(E^),  ' 


2^3  -* 


(59) 


Assuming  ° before,  we  can  equivalently  write 


(60) 


where 


A 


2-f’ 


(61) 


For  pure  capillaries  the  coefficients  are  /^2_~ 

^2  = -5o<./(32T),  while  for  pure  gravity  waves  ~3u/2  and 

0^  z 7K^c</4g^,  In  linear  theory  equilibrium  amplitudes  cannot 
exist  for  oC  < 0,  But  in  the  nonlinear  case,  a possible 
equilibrium  energy  E = exists  for 

capillaries  (E  = 16T)  but  it  does  not,  to  this  order,  for  gravity 
waves.  Eq.  (60)  can  be  integrated  for  constant  ^*s,  giving 


E(E)  = 


(62) 
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where  S is  related  to  the  initial  energy.  (For  more 
complicated  problems  in  which  the  ^*s  are  time-dependent, 
integration  in  closed  form  is  possible  using  Bemoul ' i-type 
integrals, ) 

This  result  is  reminiscent  of  the  amplitude  equation  first 

I 

given  by  Stuart  (1958)  in  nonlinear  boundary  layer  stability. 
Essentially,  if  (i)  > 0,  > 0 &n  equilibrium  energy 

exists,  (ii)  > 0,  4,  0 esi  equilibrium  energy  doesn’t 

exist,  and  (iii)  disturbance  takes  on  an 

equilibrium  energy  at  t = -«  | however,  the  equilibrium  is 

unstable.  If  the  energy  density  is  less  than  E^,  the  disturbance 
decays  to  zero,  whereas  if  it  is  greater  than  E^,  the  disturbance 
grows.  Thus,  Stuart's  equation  describes  a broader  class  of 
temporally  varying  problems.  In  other  problems  in  which  the  waves 
are  not  purely  steady  in  space  and  time,  the  modulations  balance 
according  to  Eq,  (22).  Thus  the  Landau  coefficients  in  Eq.  (62) 
are  now  functions  of  space  and  time  as  determined  from  the  wave 
kinematics.  Implicit  also  in  Eq,  (22),  of  course,  is 
Gaster's  (1962)  theorem  relating  temporal  and  spatial  growth 
rates  via  the  group  velocity!  the  relationship  is  much  more 
complicated  in  the  nonlinear  case. 


I 
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Example  6,  A Weakly  Nonlinear  Instability. 

To  examine  the  basic  effect  of  nonlinearity  on  linear 
wave  flows »•  consider  homogeneous,  conservative  media  and 
in  Eq.  (22)  with  ? = 0,  let  ua  eliminate  the  explicit  dependence 
of  and  on  A by  using  the  dispersion  relation.  This 
leads  to 

^ (63) 


which  can  be  written  as  the  following  eikonal  equation  for  the 
phase  function  5 (x,t). 


(64) 


We  can  specialize  our  discussion  to  wavemakers.  In  linear 
theory  the  wave  frequency  is  fixed  on  prescribing  the  wavenumber. 

But  in  nonlinear  theory  this  alone  is  not  sufficient  since 
something  about  the  amplitude  distribution  must  be  known.  We 
can,  of  course,  equivalently  modulate  the  frequency  and  take 
K = and  CO  s - 6 exp(lo<.x)  as  initial  conditions,  where 
£ and  o<i  are  prescribed  constants.  Then  the  perturbation  phase 
function  j>  - $ - KqX  + u/,t  satisfies  ~ ^ " £exp(io<x) 

at  t=0  as  well  as  the  differential  equation  above.  Nonlinear 
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perturbations  to  the  basic  system  thus  satisfy,  approximately, 


-li 


t± 


= o 


(65) 


The  solution  to  this  is  simply  the  real  part  of 


d ^ K^X-  U),t 


s e 


s^h  ft 
i 


(66) 


where,  if  » 


I „ C^wWjO  ijCKjO  ) 

7 


•A 


(67) 


The  instability  implied  by  Eq.  (66)  arises  from  the  solution  of 
an  elliptic  Cauchy  problem.  It  is  essentially  Whitham's  (196?) 
instability  and,  for  example,  can  be  realized  by  gravity  or 
capillary  waves  on  deep  water.  One  obvious  question  is  the 
effect  of  inhomogeniety  (e,g,,  a slowly  varying  current 
U = U(X,T)  ) and  dissipation  in  the  medium.  To  explore  this 
issue  let  US  follow  Hayes  (1970)  and  expand  the  average 
Lagrarngian  Ji  in  the  form  A - H(K,A,X,T),  This 

representation  is  a canonical  one  since  the  only  amplitude 
measure  that  appears  is  the  wave  action  A = itself.  The 
modulation  equations  as  expressed  in  Eq,  (22)  become 


I . , .1  "dK  , ij 

^kk3X  ^kx 
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(68) 


3T 


, 9K 

■^^Kax 


■))A 


(69) 


where  the  frequency  UJ  is  determineci  from 


that  is, 


w “ H^(K,A,J(,T)  (70) 

/ 

Since  is  an  energy  density  (see  Eq.  (25)),  the 

function  H can  be  interpreted  as  a Hamiltonian  energy  function. 

In  the  present  case  (as  in  Eq,  (22))  the  frequency  is  again 

reali  nonconservative  effects  appear  idrectly  in  Eq.  (68), 

As  in  the  discussion  leading  up  to  Eq,  (65)  we  consider  weakly 

nonlinear  perturbations  to  a uniform  state  (denoted  by  zero 

- subscripts)  and  linearize  Eqs,  (68)  and  (69)  about  the  constants 

and  A . If  we  introduce  new  coordinates  Z - t and 
0 o 

f = X - C-t  where  C.  is  the  "mean  velocity"  = H.„(K  ,A  »X,T), 

^00  K AK  O 0 


the  perturbation  problem  for  A*  (where  A' 
can  be  put  in  the  form 


3Z^ 


U 


2«Cu)i-Ux) 


A-Ag  and  K'  = K-K^) 


A'(i,-C-0)  - (jf-Cn 

Hci.x-o)  ® K',(i,X‘0)  « 


(71) 
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where  ( U//-  U ) is  taken  as  constant  locally.  If  £ is  email 
‘ X 

we  can  write,  approximately, 


£(ty;-Vi)X 

e 


i>) 


(72) 


In  the  hyperbolic  case  the  solution  to  the  "conservative  part" 
B satisfies 


b(i,t-o)  » (pa)  = 1/1(0 


(73) 


and  is  given  by  D'Alembert’s  formula,  that  is, 


3(\,r) 


x)+<fa-K,,KKZ«  ^ 

\ >f(V  ‘^1 


(7^) 


Eq.  (7^)  shows  how  the  perturbation  "splits" i it  also  illustrates 
the  interplay  between  • initial  conditions  and 

norj. inear ity.  The  latter  is  stabilizing  in  thesense  that  small 
perturbations  remain  small. 

The  elliptic  initial  value  problem  is  solved  by  analytic 

continuation  into  the  complex  | - plane.  The  formulation  for  B 

as  given  in  Eq.  (73)  still  holds,  but  to  indicate  explicitly 

that  fl  Htrv  0»  we  will  introduce  the  barred  variables 
A 


I 


5^ 


denoted  by 


? 


i A- J1 


a’^kk*' 


and  T = XT 


This  leads  to 


V3(I,f3  ^ U 


3(i,t=o;  “ <?)([) 


T ^ o 


(75) 


As  suggested  in  Garabedian  (1964),  we  introduce  the  complex 
variable  1 = / +-  iTf  and  assume  that  the  functions  <^  { \ ) 
and  ^ ( f ) are  analytical  functions  of  f , We  can  then  write 
B(  f , X ) = Bif-i-Cri  , T ),  B(  I , 0)  = 

and  B-(  f,T=0)  = (/'(/^tT^)  = where,  for  each  value 

of  / , B(  f , X ) satisfies 


0 


(76) 


Because  Eq,  (76)  is  formally  hyperbolic,  the  solution 
B(  X • ) can  be  determined  in  a stable  manner  from  initial 

conditions.  However,  this  is  a solution  for  complex  ? . The 
physical  solution  is  obtained  in  the  limit  0 after  an 

explicit  formula  for  B(  i , T ) is  obtained.  Again,  the  D'Alembert 
formula  applies,  giving 


(77) 
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where  the  semicolon  above  refers  to  the  parametric  dependence 
• oa  ^ (the  physical  solution  implies  the  limit  I -S  ).  As  a 
simple  example  consider  the  situation  where  the  initial 
wavenumber  variations  are  zero  so  that  7/  ) = 0,  Since 

Eq.  (76)  is  linear  it  suffices  to  consider  a Fourier  component 
B(T=0,  I ) = oL  ain^t  of  the  initial  action  disturbance,  that 
is,  ■ oC  ain  ^ , In  this  case,  the  physical 

solution  is 

?l-*o  ^ 

The  appearance  of  the  cosh^^  in  addition  to  the  growth  rates 
U)i  and  already  in  Eq,  (72)  enable  a simple  definition  for  an 
"effective  c<y^.". 

In  applying  D'Alembert's  formula  and  Garabedian's  method 
of  imaginary  characteristics  certain  conditions  on  smoothness  and 
analyticity,  of  course,  are  implicitly  assumed.  These  are 
in  addition  to  those  calling  for  slow  variation,  as  indicated  in 


Eq,  (2),  Sufficiently  non-smooth  Cauchy  data  also  possess  real 
solutions,  but  solutions  to  these  problems  are  not  likely  to 
be  slowly  varying,  and  hence,  are  not  likely  to  be  described  by 
the  present  theory. 
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III,  High-Order  Effects. 

The  low-order  formulation  discussed  in  Section  II 
gives  results  which,  strictly  speaking,  are  invalid  over 
space/time  scales  greater  than  0(6'*),  They  are,  of  course, 
valid  in  the  sense  under  which  they  were  derived,  but  to  deal 
more  accurately  with  large-scale  effects,  high-order  modulations 
as  described  in  Eqs,  (1?)  and  (20)  must  be  considered, 

Lighthill  (1965a),  as  an  example,  traced  the  progress  of  gravity 
wave  packets  on  deep  water  and  showed  how,  on  the  basis  of 
low-order  Whitham  theory,  cusp-like  energy  redistributions 
within  the  packets  would  asymptotically  result.  However,  the 
inclusion  of  high-order  dispersive  terms  (Yuen  and  Lake,  1975) 
actually  suppresses  the  Lighthill  singularity!  the  end  result 
is  a train  of  solitons.  Another  situation  where  high-order 
dispersive  and  diffusive  corrections  are  likely  to  be  important 
is  in  the  formation  of  so-called  "kinematic  shocks".  These  are 
essentially  wave-mechanical  discontinuities  which,  on  the  basis 
of  low-order  theory,  are  assumed  to  conserve  various  fluxes, 
for  example,  momentum  or  energy.  Their  probable  existence  was 
first  speculated  upon  by  Whitham  (I965),  However,  the  true 
nature  and  stability  of  these  shocks,  of  course,  cannot  be 
studied  without  considering  the  structure  of  the  high-order  terms. 
It  is  entirely  possible  that  sufficiently  strong  dispersive 
effects  in  the  high-order  correction  terms  can  suppress  their 
ever  forming.  On  the  other  hand  if  the  corrections  are 
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essentially  diffusive,  as  in  gas-dynamic  shocks,  the  low-order 
results  would  be  meaningful  in  the  sense  that  the  solution  is 
"properly  embedded"  within  the  high-order  one.  More  on  this 
will  be  said  in  Section  IV, 

For  many  applications  in  continuum  physics  the  basic  wave 
solution  is  weakly  nonlinear,  and  hence,  expandable  via  Fourier 
series  in  powers  of  amplitude  squared.  For  these  cases  the 
structural  form  of  the  high-order  modulation  terms  depends 
only  on  the  primary  harmonic,  that  is,  the  basic  linear  wave 
solution.  In  this  physical  limit  one  can  obtain  a more 
transparent  but  completely  equivalent  re-expression  of 
Eqs,  (17)  and  (20),  in  which  the  corrections  to  low-order  theory 
depend  only  on  the  real  and  imaginary  frequencies  jZ(K)  and 
J^‘'(K)  corresponding  to  the  uniform  wave  solution.  To 
facilitate  the  discussion,  let  us  first  consider  purely  linear 
problems  without  flow  inhomogenieties , so  that  Fourier  super- 
position is  valid.  The  mathematical  approach  is  a simple  one. 
Expansions  about  a center  wavenumber  and  a real  frequency  Cd^ 
for  the  modulations  are  introduced  as,  for  example,  in 
Davey  (1972),  but  unlike  Davey  the  series  solution  is  not. 
truncated  to  a finite  number  of  terms.  The  full  series  solution 
is  then  resummed  in  a type  of  analytic  continuation  to  produce 
solutions  valid  for  all  wavenumbers  that  are  consistent  with  the 

Lagrangian  approach.  The  technique  used  is  analogous  to,  for 

•o 

example,  identifying  the  series  ^ with  (1-x)  , which 

nso 

has  a wider  range  of  convergence.  We  now  consider  the  linear 
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superposition  of  monochromatic  wave  components  exp  iiKx-  UA) , 
each  satisfying  the  complex  dispersion  relation 


UJ  = J1(K)  = Jl" 


*iJl‘ 


(79) 


Here  K is  a real  wavenumber  with  Jl.^  and  _fl  being  respectively 
the  real  and  imaginary  parts  of  „/l(K),  In  the  neighborhood  of 
the  center  wavenumber  K^,  Taylor  expansion  gives  the  series 

2 . If  ^ 0(  & ) « 1 we  can,  as 

in  Lighthill  (1965b),  Fourier  superpose  elementary  solutions 
and  arrive  at  the  more  general  solution 

/V.  r”  UKX-u/i) 

F(X,i)  - J F(K)  e (80) 


The  function  B(K)  is  proportional  to  the  Fourier  transform  of 
the  initial  condition  and  it  exists  when  the  initial  disturbances 
are  localized.  Without  loss  of  generality  let  us  re-express 
Eq,  (80)  in  the  form 


F(x,i)  - ^(x,i)  e 


e 


(81) 


(82) 
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80  that  F(x,t)  consists  of  a purely  periodic  part  and  an 
amplitude  function  ^(x,t)  that  includes  dissipative  effects. 
This  follows  the  approach  introduced  in  Section  II.  Applying 

the  operator  J ...  0 d K to  the 

— 

Taylor-expanded  dispersion  relation  (which  holds  for  any 
component  wave)  and  noting  that  L ^ 

and  (ri)  Px"  * Q- leads  to  the  following 

identity  for  . 


(83) 


Next,  as  in  Eq.  (2),  we  seek  WKB  solutions  of  the  form 


f * aiX.V  e 


Le(x,t) 


(84) 


where  ^ is  the  perturbation  phase  function,  that  is, 

=X(X.T)  and  Zl7(X,T),  Now  the  multiple- 

scaling  introduced  in  Eq.  (2)  implies  that,  for  any  function 
u(x,t)  = U(  (5  ,X,T), 
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u„{x,t; 


I"U„5 

+£^  . nfr)-iKM'2)('1-3;  yrf^-A-  ] y 

Q h.  ^^^(n-2)e 

. yifn-0(»-^)  IT  V 1 r 

g li.  ^(n-2)9 

ii(M-i)(M-z)(n-:3)fi-+}  ^»-5^  y ,_ 

/2.  Ji  ^x-^XX 

. nfn~i)(^-Z)(in-i)  ir^-^rr  7r 

+ g A X 

. Y\(lr)~0(^~V(^'^^  IT  nr 

3 A i^x^(t^-})s,y.X 

. n(n-iXn-2)  ■i^''-^  ? r 

6 A.  ^n-3)9,XXX 

J.  trj(i1-0(H'2X>l'3-?  -cr  T7-  7r 

2?  ^ 

, nCy)’\)(*\-z)[rt-3Xy\'^Xri-5)  ^n-4  3 ,, 

^ 43 L Xx  ^Ln-3)d 


(85) 
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aa  can  be  readily  shovm  by  induction.  Let  us  expand  the 
function  (9»X,T)  in  Eq,  (84)  following  Eq,  (85)  and  substitute 
the  results  into  Eq.  (83),  remembering  that 
satisfiec  1.  One  sees  that  to  obtain  0(2’’ ) 

corrections  to  the  linear  phase  relation  only  0(£.’’)  terms  in 
Eq.  (85)  need  to  be  retained.  But  to  obtain  O(t^)  corrections 
to  the  low-order  amplitude  equation  (modified  to  include  the 
effects  of  dissipation)  we  must  retain  the  0(£^)  terms  in 
Eq.  (85),  Equating  real  and  imaginary  parts  in  the  equation 
just  obtained  gives,  respectively, 
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u‘jtA-4.. 


-^iiK.o 


n*i  Hi 
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il»l  Hi 


-eV 


K'Xx  " ^ 


:x  X 


+ 6^1'!^® 
n*i  n I 


'1C>1-I)(H-2)(I7-3XI1-+^  T7-"-*  T7-  ir 

ii  JS^  ^x^xx"- 
G J1  ^xx^x  Sxx 

4.  rtf»-/X»?-2)(H'3;-prn-+p-  ^ 

+ ^ ^X^xx 

^ l1fHH}CH-2)(l7-3;  i^n-^  ^ 

Zf  ^ L 

^ w(ii-iXn-2)(n-3Xw-4-;  ^ 


X X 


+ )(H*'3^C^~4<^(I'?*  5/)  ^ W"  6 ^ 3 ^ 

^8  X 


The  first  equation  describes  the  perturbation  phase 
corresponding  to  the  slowly  varying  wavetrain  and  the  second  one 
describes  amplitude.  However,  since  we  expect  our  results  to 
extend  to  all  wavenumbers,  we  sum  the  above  series.  This  leads 
to  a phase  relation  of  the  form 
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60=il'-^e-ai  ^ 


a 


I^K  2 


?2 

'KK  X 


(88) 


After  some  manipulation  the  amplitude  equation,  in  X-T 
variables,  becomes 


= 2f 


(89) 


-e* 
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Si 
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KK 


+ aa^K^  f 


XX  £ 


^^XXX-^3K  '^2  ^x'^*XX-^4<r 
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Because  il  //I  0(  E ) all  high-order  corrections  to  the 

basic  low-order  phase  are  0(£^),  as  would  be  expected  for 
weakly  nonconservative  waves  that  are  truly  slowly  varying. 

(This  is  also  consistent  with  the  approximate-Pourier  superposition 
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used.)  In  the  Lagrangian  approach,  the  vanishing  of  0( S ) 
terras  would  be  attributable  to  the  orthogonality  of  the 
trigonoraetric  functions,  and  identical  results  (Chin,  1976) 
can  be  obtained  by  extending  the  method  of  Example  1 in  the 
manner  of  Yuen  and  Lake  (1975)i  the  latter  of  which  will  be  later 
discussed.  Similarly,  all  high-order  corrections  to  the 
low-order  amplitude  equation  with  damping  are  0(£  ).  In  general, 
the  indicated  corrections  involve  coefficients  that  depend  only 
on  high-order  wavenumber  derivatives  of  Sif'  and  -TL  , which  are 
in  turn  defined  by  the  planar  monochromatic  wave  solution  for  a 
particular  problem. 

The  above  results  show  how,  in  a simple  way,  the  complex 
dispersion  relation  of  the  monochromatic  wave  completely 
determines  the  dynamics  of  the  slowly  varying  one,  as  would  be 
expected  from  the  physics  of  the  problem.  The  high-order  terms 
in  Eqs.  (88)  and  (89)  based  on  Jl.  (K)  and  its  wavenumber 

derivatives  are,  of  course,  directly  related  to  the  high-order 

« 

Lagrangian  terms  in  Eqs.  (17)  and  (20)j  those  based  on  Jl. 
appear  on  account  of  the  dissipation-,  functional  F,  But  the 
Fourier  integral  approach  given  here  provides  a much  more 
transparent  indication  of  the  structure  of  the  high-order 
dispersive  and  diffusive  effects  than  does  the  equivalent 
variational  onej  the  results  derived  here  are  much  simpler  to 
use  and  do  not  bear  the  limiting  assumptions  central  to  other 
approaches  employing  centered  wavenumber  expansions.  Let  us 
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now  consider  the  effect  of  weak  nonlinearity.  As  previously 
discussed  the  high-order  modulation  terras  for  a weakly  nonlinear 
system  are  determined  from  the  primary  harmonic,  that  is,  they 
are  as  given  in  Eqs.  (88)  and  (89).  Thus,  nonlinear  Stokes-type 
corrections  can  be  incorporated  in  the  above  results  in  a 
simple  "additive"  manner.  For  simplicity  let  us  first  deal 
with  conservative  problems.  The  low-order  amplitude  equation 
is  just 


(90) 


where  U)  =.Jl  (K)  is  the  linear  dispersion  relation.  For 
moderately  small  amplitudes  a Stokes  expansion  implies  that 


CO  ==  jta)  + 4co  d'' 


(91) 


and  the  wave  conservation  condition  dK/dt  + 9 w//^x  = 0 can 

be  expanded  out  to  give 


(92) 
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The  important  coupling  arises  in  f2(K)  8a  /8x  since  it  implies 
an  0(a)  correction  in  the  characteristic  velocities.  The  other 
new  term  merely  corrects  the  coefficient  of  the  existing  term 
in  3K/^x  and  contributes  only  at  the  0(a  ) level.  Similarly, 

the  nonlinear  corrections  to  Eq.  (90)  would  be  various  terms 

4 2 

of  0(a  ) providing  relative  corrections  of  0(a  ) to  the 

coefficients  of  the  existing  terras  in  8a  /8x  and  8K/ 9x. 

Thus,  in  the  first  assessment  of  nonlinear  effects,  we  can 

retain  Eq.  (90)  and  use  the  new  dispersion  relation  as  given 

in  Eq.  (91),  dropping  the  nonlinear  correction  to  ^ in 

Eq.  (92).  This  is,  in  essence,  Whitham's  (I967)  approximation. 

But  since  high-order  modulation  effects  arising  from  linear 

terms  must  eventually  interact  with  the  nonlinearity,  Whitham’s 

approximation  should  be  modified  to  describe  this  effect. 

As  is  clear  from  the  above  discussion,  the  leading  high-order 

Corrections  will  modify  only  the  phase.  Bearing  in  mind  that 

for  near-linear  problems  the  corrections  are  "additive"  in  the 

sense  previously  discussed  we  have,  on  assuming  that  wavenumber 

variations  are  less  rapid  than  those  in  amplitude  (or,  if 

3K  • '^4'<  » •••  neglected)  the  following  approximation 

w - i2’'cK)  + -f.ma’-  - i ^ (93) 

whd^h  must  be  solved  together  with  Eq.  (90),  These  equations, 
again,  hold  only  for  homogeneous,  conservative  media. 
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To  understand  the  qualitative  nature  of  the  above 
approximation,  let  us  re-expand  Eqs,  (90)  and  (93)  about 
K = Kq  , say,  consistently  with  our  assumption  of  slow  wavenumber 
variations.  Again  considering  the  conservative  case  we  have, 
with  Co  = Jl'  o ■ 

+ (Tc.+ii'ic^,  =■  o (91*) 

U3-(M,  =C.(<-K.)  (95) 

<v 

Using  the  notation  ^ = ae^*^  , these  equations  combine  to  give 

“/zCO /)('/>  (96) 

This  "nonlinear  Schroedinger  equation"  has  been  derived  in  a 
number  of  contexts  by  various  authors  using  different  techniques. 
But  our  discussion  shows  that  it  is  really  nothing  more  than  a 
specialized  form  of  Eqs.  (1?)  and  (20).  This,  in  fact,  will 
be  clearly  illustrated  in  Example  7. 

Now  in  "self -focusing"  problems,  Zakharov  and  Shabat  (1972) 
have  shown  that  for  initial  conditions  which  approach  zero 
sufficiently  rapidly  as  1x1""®®  (corresponding  to  "pulses"), 

Eq.  (96)  can  be  solved  exactly  by  the  inverse  scattering 
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I technique.  They  discovered  that  an  initial  wave  envelope 

pulse  of  arbitrary  shape  will  eventually  disintegrate  into  a 

number  of  solitons  and  an  unsteady  oscillatory  "tail"  interpreted 

as  noise.  The  number  and  structure  of  these  solitons  and  the 

structure  of  the  tail  are  completely  determined  by  the  initial 

conditions.  But  the  tail  is  relatively  small  and  unimportant 

for  pulse-type  initial  conditions;  it  disperses  linearly, 

resulting  in  a amplitude  decay.  Each  of  the  resulting 

solitons  is  a permanent  progressive  wave  characterized  by 

(four)  parameters  in  amplitude,  speed,  position  and  phase,  but 

unlike  the  solitary  wave  solutions  of  the  KdV  equation  for 

shallow  water,  the  amplitude  and  speed  parameters  bear  no 

direct  relation  to  each  other  except  that  they  are  the  imaginary 

and  real  parts  of  an  associated  scattering  problem.  It  is 

important,  though,  that  these  solitons  are  stable  in  the  sense 

that  they  can  survive  interactions  with  each  other  without 

permanent  change,  except  a possible  shift  in  (the  parameters 

for)  position  and  phase.  The  time  scale  for  the  soliton 

formation  was  found  to  be  in  direct  proportion  to  the  length  of 

the  pulse  and  in  inverse  proportion  to  the  amplitude  of  the  pulse. 

The  above  result  shows  that  the  end  product  of  unstable 

"elliptic"  modulations,  unlike  Lighthill's  (1965a)  low-order 

cusp  solution,  is  a train  of  solitons.  The  exact  behaviour 

depends  on  the  relative  magnitudes  of  and  f«(K^),  but 

IVA  ,0  c 0 

the  same  qualitative  features  can  be  expected  for  unstable 
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gravity  and  capillary  waves,  or  for  that  matter,  all  waves 
that  are  "elliptically  unstable”  in  Whitham's  (196?)  sense.  To 
show  that  the  approximation  underlying  Eqs.  (90)  and  (93).  in 
fact,  can  be  obtained  by  Lagrangian  techniques,  it  suffices  for 
our  purposes  to  present  one  illustration. 


Example  7,  Gravity  Waves  in  Deep  Water, 

We  quote  the  example  considered  by  Yuen  and  Lake  (1975). 
who  chose  for  the  basic  Ansatz  (corresponding  to  Eq.  (21))  the 
wave  solution,  correct  to  0(  £ ,Ka), 

71  aCEX.ti)  cos  6 CoslB  * ...  (97) 

^ ^ ^ cos  5 (l-tC^JcoseJc  sin  18*^  ( 98 ) 


The  Lagrangian  density  is  just  the  pressure,  that  is. 


where  ^(x,t)  is  the  free-surface  position,  y is  the  modal 
coordinate,  <p  (x,y,t)  is  the  velocity  potential,  g is  the 
acceleration  due  to  gravity,  with  the  fluid  density  being  unity. 
Then,  the  average  Lagrangian,  correct  to  0( £ ,K  a ),  is 
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^ -i- 

‘ ^K  4-  ^K  4-K^  ZK  (100) 

z ^ 

^ ^ ^ _3  UJ^U^XX  ^3^ 

4-  K^  B i^S  8 


and  we  may  note  that,  in  this  approximation,  space-time 
derivatives  of  frequency  and  wavenumber  do  not  appear.  After 
some  algebra,  they  found  that  S variations  ga-'e 


H 


PA  *"3 


O 


(101) 


while  "a"  variations  gave,  on  simplifying,  the  phase  relation 


UJ 


= H- 


^xx  -I 

8clK^  J 


(102) 


Now  this  is  precisely  the  approximation  deduced  in 
Eqs.  (90)  and  (93).  To  see  this,  simply  note  that  the  Stokes 
solution  gives  f2(K)  = i g^  while  -^ysl  ~ 

from  the  linear  part  of  tha  dispersion  relation.  The  above 
results  were  also  derived  by  Chu  and  Mei  (1970)  using  multiple- 
scaling  methods  directly  from  the  governing  equations.  All 
these  results,  of  course,  are  really  special  cases  of 
Eqs,  (88)  and  (89)  modified  for  weak  nonlinearity,  which  do  not 
bear  the  limitations  imposed  by  "centered  wavenumber"  expansions. 
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This  is  especially  important  for,  in  analyses  of  the  latter 
kind,  one  implicitly  assumes  (in  some  nondimens ional  sense)  that 
amplitude  variations  dominate  those  in  wavenumber  - otherwise 
the  expansions  are  meaningless.  In  another  sense,  the 
generalized  modulation  equations  derived  are  useful  because 
it  is  not  even  clear  that  a "centered  wavenumber"  stays 
constant  on  the  scale  over  which  the  high-order  modulation 
terms  become  important.  It  is  crucial  to  note  the  equivalence 
between  Sqs.  (20)  and  (1?)  and  Eqs.  (88)  and  (89).  They 
remain  applicable  for  large  values  of  K^(X,T). 

To  the  above  results  for  conservative  and  homogeneous 

media  we  now  add  the  effects  of  dissipation.  The  same  line  of 

reasoning  that  led  to  Eq.  (93)  In  the  nondissipative  case  now 

requires  that  we  retain  in  Eq.  (88)  the  term,  which 

balances  the  dispersion  term,  both  being  0(£^).  At  the  same  ' 

time  the  low-order  damping  term  2'^  in  Eq.  (89),  first 

suggested  in  Landahl  (1972),  must  be  kept  (see  Example  1). 

The  first  of  these  corrections  seems  to  have  appeared  in  a 

number  of  studies  in  hydromagnet ic  waves,  according  to 

Malkus  (1976),  though  in  a less  general  form.  In  the  next 

approximation,  high-order  corrections  to  the  "simply-damped" 

Ji!" 

amplitude  equation  involve  the  diffusion  terras , 
etc,,  and  the  dispersion  terms  » etc.,  in  Eq.  (89), 

and  the  remaining  terms  in  Eq.  (88).  The  etc.’s  here  generally 
refer  to  terms  involving  higher  differentiations  of  fl  and  Si 
with  respect  to  wavenurabei'  and  the  appearance  of  the  terms 
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' ^XX  ^XXX  would  have  been  omitted  in  the  centered 

..wavenumber  approach).  Thus  the  appearance  of  these  terms 
depends,  in  a sense,  on  "how  dispersive"  or  "how  dissipative" 
the  original  uniform  wavetrain  is.  For  example,  the  vanishing 
of  Jl^l^  would  leave  only  diffusive  corrections  to  Eq.  (89) j the 
phase  relation,  though,  would  still  bear  the  effects  of 
high-order  dispersion.  In  this  kind  of  limit,  the  most  obvious 
effect  of  diffusion  would  appear  to  be  a gradual  decay  of  the 
solitons  as  are  produced  in  the  Zakharov  and  Shabat  (1972) 
solution.  No  general  conclusions  are  available,  at  this  point, 
for  Eqs,  (88)  and  (89)  with  (or  without)  their  nonlinear 
corrections.  Some  future  studies  are  planned,  however,  to  test 
the  separate  effects  of  strong  dispersion  and  strong  diffusion 
for  various  kinds  of  waves.  However  some  speculation  is  possible 
on  what  types  of  "strange"  behaviour  might  be  expected.  For 
example  the  "diffusion  coefficient"  in  Eq,  (89)  in  the 

case  of  deep-water  surface  waves  is  just  the  kinematic 
viscosity  (Lamb,  1932)  which  is  positive,  as  expected.  This 
positivity  is  also  apparent  from  existing  numerical  studies  in 
boundary  layer  and  vortex  instability  1 to  the  author's 
knowledge,  the  plot  of  J)! {K)  against  K has  always  appeared 
concave  down.  A surprising  result,  however,  is  found  in  warm 
electron  plasma  waves.  Classically  the  waves  are  Landau-damped, 
that  is,  -Ti  (K)  ^ 0,  so  that  according  to  the  low-order  form 
of  Eq.  (89),  the  waves  always  decay.  However,  direct  calculation 
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of  shows  that  there  exists  a range  of  wavenumbers 

bearing  negative  diffusion.  Thus,  on  a large  time  scale,  there 

may  occur  a transfer  of  energy  back  into  certain  of  the  formerly 

decayed  waves.  The  possibility  of  sign  variations  in  - _/2^ 

KK 

is  of  fundamental  importance,  especially  from  the  viewpoint  of 
developing  deterministic  wave  models  of  turbulent  shear 
flow  (Chin,  1976)}  it  allows,  for  example,  energy  transfer  by 
classical  cascading  and  by  reverse-cascading. 

We  also  mention,  in  context,  the  "nonlinear  instability 
burst"  found  by  Stewartson  and  Stuart  (197I),  who  studied  the 
stability  of  wave  systems  in  plane  Poiseuille  flow.  They 
considered  an  infinitesimal  centered  disturbance  imposed  on  a 
fully  developed  plane  Poiseuille  flow  at  a Reynolds  number  R 
slightly  greater  than  the  critical  value  R^  for  instability. 

The  disturbance  was  assumed  as  a wave  modulated  in  space  and 
time,  whwse  amplitude  "a"  satisfied  a nonlinear  parabolic 
equation.  For  finite  values  of  (R-Rjj)t,  Hocking,  Stewartson 
and  Stuart  (1972)  showed  that  the  amplitude  develops  an  infinite 
peak  at  the  center  of  the  wave  group.  This  result  also  appears 
an  a special  limit  of  the  present  theory,  although  the  number 
of  assumptions  required  seems  somewhat  restrictive.  In 
Stewartson* s approach,  center  wavenumber  expansions  for  K = 
are  taken  about  a point  of  "maximum  growth",  that  is, 

= 0,  If  we  further  assume  that  the  waves  are  weakly 
dispersive  - 0)  the  phase  relation  becomes  = J2.  (k) 

and  the  amplitude  equation  reduces  to 
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2T  ^ t 


(103) 


Stewartson  then  determined  the  coefficients  of  "a"  and  its 
derivatives  on  the  basis  of  linear  stability  theory  and  added 
O(a^)  nonlinear  corrections  to  the  growth  rate  jQ.^(K)  by 
expanding  the  production  terra  in  powers  of  amplitude.  This  led 
to  a nonlinear  parabolic  equation,  whose  solution  was  obtained 
by  matching  to  the  stationary  phase  solution  of  linear 
theory  (assuming  zero  disturbances  far  from  the  wave  center). 

The  present  author's  own  view  is  that  Stewartson' s theory 
(as  generalized  here)  is  overly  restrictive  in  the  sense  that 
the  kinematic  coupling  between  wavenumber  and  wave  amplitude  is 
not  considered,  This  may  be  correct  in  the  initial  stage,  but 
it  is  certainly  not  so  in  the  latter  stages  during  which  the 
amplitude  singularity  appears.  We  mention  in  passing  that  the 
structure  of  Landahl's  (1972)  focus  (based  on  the  effects  of 
ray  focusing  due  to  inhomogenieties  using  Whitham’s  equations 
modified  for  dissipation)  was  also  obtained  in  Chin  (1976). 

Matched  asymptotic  expansions  were  used  to  connect  the  "inner" 
focal  structure  (based  on  a modified  form  of  Eq.  (89)  allowing 
for  inhomogeniety ) where  diffusive  effects  would  be  important 
to  an  "outer"  solution  based  on  Whithara's  equations  modified 
to  account  for  nonconservative  effects.  The  waves  were  assumed 
(1)  to  be  only  weakly  dispersive,  (2)  to  have  ~ approximately 

constant,  and  (3)  to  satisfy  (dA/A )/(dK/K) » 1 near  the  focus. 
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As  was  expected,  the  effect  of  diffusion  was  a smoothing  of 
the  amplitude  peak.  The  thickness  of  the  region  over  which 
diffusive  effects  are  likely  to  be  importnat  was  found  to 
be  of  ) in  unsealed  x-space.  This  thickness  decreases 

with  decreasing  viscosity  8 , as  would  be  expected  on  physical 
grounds.  In  contrast,  Stewartson's  "nonlinear  instability 
burst"  for  channel  flows  occupies  a length  of  ) 

which  increases  with  decreasing  viscosity.  The  details  of  the 
above  study  will  be  presented  in  a separate  paper. 
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IV.  Discontinuous  Solutions. 

^ Perhaps  the  most  fascinating  area  of  nonlinear  wave 
mechanics  is  the  study  of  shocks  and  their  formation.  The 
possible  existence  of  these  wave  mechanical  discontinuities 
was  first  suggested  by  Whitham  (1965)  and  a good  discussion 
on  the  subject  appears  in  Chapter  15  of  Whitham  (197^). 

However  there  is  little  experimental  evidence  as  to  their  actual 
existence.  The  observability  of  the  shocks  suggested  by 
Whitham* 3 low-order  theory  may,  in  practice,  be  somewhat 
obscured  by  the  effect  of  high-order  dispersive  and  diffusive 
terms  as  were  previously  derived  in  Section  III  and  the  effects 
of  low-order  inhoraogeniety,  which  are  certain  to  affect  their 
stability.  The  first  type  of  shock  considered  by  Whitham  deals 
with  the  question  of  "breaking"  and  arises  when  the  low-order 
modulation  equations  are  hyperbolic.  The  dependence  of  the 
characteristic  velocities  on  the  modulation  variables  introduces 
the  usual  hyperbolic  distortion,  and  "compressive"  modulations 
in  a simple  wave  solution  will  develop  multivalued  regions. 

It  is  possible,  Whitham  notes,  that  these  solutions  actually 
represent  superpositions  of  two  or  more  wavetrains  with  different 
ranges  of  K and  "a".  The  solutions,  though,  would  not  be 
correctly  described  by  Eq.  (2)i  some  type  of  multiphase 
generalization  would  be  needed,  Whitham  compares  this  situation 
to  that  in  linear  theory  where  the  group  velocity  C_(K)  decreases 
toward  the  front.  Then,  since  values  of  K propagate  with 
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velocity  C (K),  some  kind  of  overlapping  would  occur.  In  the 

O 

linear  case  the  complete  process  is  described  by  Fourier 
integrals?  presumably,  in  the  nonlinear  case,  some  type  of 
nonlinear  superposition  would  hold.  There  does,  however, 
exist  the  possibility  that  higher  order  terras  become  important 
near  breaking  and  prevent  the  development  of  multivalued 
solutions.  This  is  certainly  the  case  if  the  only  ) 
corrections  kept  are  the  dispersive  ones,  that  is,  those 
dependent  on  (K)  as  shown  in  Eqs.  (88)  and  (89).  The 
resulting  equations  then  become  similar  in  form  to  the  Boussinesq 
and  Korteweg-deVries  equations,  and  by  analogy,  breaking  would 
be  suppressed.  This  behaviour  would  most  likely  apply  to 
small  symmetric  modulations,  as  noted  by  Whitham  (197^).  and 
eventual  development  into  series  of  solitary  waves  would  occur, 
as  would  be  suggested  by  the  solution  of  Zakharov  and  Shabat  (1972). 
On  the  other  hand,  strong  unsymmetrical  modulations  may  break 
in  some  sense.  If  the  high-order  diffusion  terms  indicated  in 
Eqs.  (88)  and  (89)  are  more  dominant  than  the  dispersive  ones, 
the  resulting  amplitude  equation  resembles  a modified  Burgers* 
equation.  In  the  narrow  "breaking  region"  where  diffusion 
balances  inertial  effects,  the  J2  expected 

to  produce  a continuous  shock  structure.  It  is  interesting  to 
note  that  relative  effects  between  high-order  diffusion  and 
dispersion  will  in  general  be  different  for  different  wavenumbers. 
Thus,  the  situation  will  be  much  more  complicated  than  the 
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analogous  case  of  gasdynamic  shocks.  There  is  still  another 
important  difference.  Provided  - is  nonzero,  it  is 

always  possible  to  define  a diffusion  coefficient  even  for 

a perfectly  inviscid  system.  For  surface  waves  on  deep  water 
there  is  no  difficulty!  the  diffusion  coefficient  so  defined 
is  just  the  kinematic  viscosity.  But  for  waves  on  inviscid 
shear  flows,  wave  diffusion  will  be  present  in  general.  It  is 
even  possible  for  this  diffusion  to  be  negative  in  a certain 
range  of  wavenumbers  (as  was  the  case  for  waves  on  warm  electron 
plasmas).  Then  the  gasdynamic  analogy  calling  for  a smooth 
shock  structure  would  not  hold.  It  is  not  clear,  though,  what 
would  happen;  but  the  high-order  equations  derived  in  Section  III 
(suitably  modified  for  nonlinearity)  can,  in  principle,  be 
considered  numerically  as  initial  value  problems. 

The  second  class  of  shocks  posed  by  Whitham  arises  in  the 
search  for  weak  solutions  to  the  modulation  equations.  For 
problems  in  which  a mean  background  state  does  not  enter 
(e»g.»  gravity  waves  in  deep  water),  the  only  dependent  variables 
are  those  for  the  wavenumber  K and  the  wave  amplitude  "a". 

The  values  of  "a"  and  K upstream  and  downstream  of  the  discontinuity 
would  be  connected  by  postulated  jump  conditions,  for  example, 
those  conserving  energy  flux,  momentum  flux,  frequency,  etc. 

The  actual  choice  of  jump  conditions  is  open  to  question, 
however.  One  can  argue  that  energy  and  momentum  are  conserved  in 
the  detailed  description  for  the  variable  u (see  Eqs.  (1)  and 
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(2))  and  should  therefore  be  retained  in  the  slowly  varying 

approximation.  These  shocks  would  therefore  represent  a source 

of  oscillations  and  involve  jumps  in  the  adiabatic  invariants, 

"If,  on  the  other  hand,"  quoting  Whitham  (197^),  "the 

discontinuities  are  supposed  to  represent  phenomena  not 

covered  by  the  original  equation,  but  covered  by  some  even 

more  detailed  description  involving  dissipation  of  some  kind, 

then  the  choice  would  be  different.  Although  momentum  would 

probably  be  conserved,  energy  presumably  would  not,"  Action 

flux  is  not  likely  to  be  conserved,  "but  one  could  make  a case 

for  frequency.  With  dissipation  smooth  oscillatory  changes 

between  different  constant  states  may  be  constructed  in 

dispersive  models,"  The  constant  end  states  referred  to  above 

would  be  the  result  of  dissipation  dampening  out  the  oscillations 

on  the  two  sides  of  the  transition  region.  The  qualitative 

effects  of  diffusion  were  considered  in  Whitham  (197^,  p.  ^82) 

using  a simple  model.  There  he  added  a Burgers  type  terra  to 

the  KdV  equation  to  produce  a simple  model  for  the  structure  of 

bores.  The  actual  terms,  of  course,  are  given  in  Eqs,  (88) 

and  (89)  if  the  system  is  weakly  nonlinear.  The  dissipation 

terra  is  really  ~ jT}!’  a , suggesting  that  the  unsteady 

KK 

formation  and  resulting  shock  structure  (not  attempted  here) 
is  much  more  complex  that  Whitham' s model  would  suggest, 
although  his  results  are  qualitatively  correct.  The  actual 
results,  of  course,  would  be  highly  wavenumber  dependent.  The 
existence  of  the  above  shocks,  at  this  point,  is  speculative. 
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But  they  would  correspond  to  the  low-order  weak  solutions 
of  the  full  equations  (in  the  appropriate  conservation 
form)  with  the  dissipation  functional  F set  identically  to  ze: 
at  the  outset.  The  high-order  terms  in  Eqs.  (88)  and  (89) 
do  provide  the  answers  as  to  the  structure  of  the  transition 
region  ^ the  above  shocks  indeed  exist.  They  also  lay  the 
groundwork  for  a first-cut  phase  plane  analysis. 
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A strong  reason,  perhaps,  for  studying  discontinuous 
solutions  of  the  above  kind  is  suggested  by  an  analogy  with 
gasdynaraics.  Consider  the  propagation  of  one-dimensional 
unsteady  disturbances  in  a transonic  decelerating  channel  flow. 

It  is  possible  to  trace  the  evolution  of  expansion  and  compression 
pulses  originating  downstream  of  a linearly  decelerating  sonic 
region  by  examining  the  equations  for  one-dimensional  unsteady 
isentropic  flow.  Kantrowitz  (1947)  was  first  to  show  that  such 
disturbances  are  deformed  as  they  propagate  upstream,  and  that 
the  disturbances  tended  to  collect  in  the  sonic  region.  The 
deformation  of  the  pulse  shape  actually  causes  weak  shock  waves 
to  form,  either  at  the  head  front  or  the  rear  of  the  pulse; 
expansion  pulses,  which  terminate  in  shock  waves  at  their  rear, 
are  subsequently  "consumed’’  by  their  own  shock  waves,  while 
compression  pulses,  which  have  shock  fronts,  grow  in  strength 
as  their  shock  fronts  progress  into  the  supersonic  region. 
Ultimately,  in  a real  channel,  they  would  destroy  the  supersonic 
flow,  converting  it  into  subsonic  flow.  This  simple  example 
suggests  the  existence  of  flows  which  cannot  support  arbitrary 
infinitesimal  perturbations,  but  which  must  break  down  locally 
into  finite  amplitude  initially  small-scale  oscillations. 
Instabilities  of  the  above  kind  are  "catastrophic"  in  nature, 
but  they  are  not  restricted  to  transonic  flows.  One  can 
speculate  that  similar  highly  unsteady,  nonlinear  physical 
mechanisms  might  be  responsible  for  other  "strong" 
fluid-dynamical  instabilities,  for  example,  as  in  hydraulic 
jump  formation,  vortex  breakdown,  and  perhaps  the  sudden 
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transitions  observed  in  boundary  layer  instability.  The 
transonic  analogy  leads  us  to  envision  an  initially  (essentially 
-linear)  phase  during  which  selected  self-excited  secondary 
waves  focus  on  a larger  scale  inhomogeneous  primary  wavetrain. 
Subsequent  trapping  and  continued  wave  growth  in  the  absence 
of  dissipative  effects  would  then  precipitate  a rapid,  almost 
discontinuous,  and  certainly  a strongly  nonlinear  readjustment 
of  the  mean  flow  in  order  to  satisfy  certain  global  constraints 
(for  example,  conservation  of  total  mass,  momentum  or  energy). 
This  type  of  breakdown  should  be  a distinct  possibility  for  any 
inhomogeneous  continuum  system  which  can  support  propagating 
waves  of  scale  small  relative  to  the  scale  of  the  inhoraogeneit ies 
The  "initial  x^Iiase"  (as  described  above)  is  amenable  to  general 
linear  analysis  using  recently  developed  ideas  in  kinematic 
wave  theory,  and  first  results  pursued  along  these  lines  have 
been  given  by  Landahl  (1972).  (Kinematic  wave  theory  is 
developed  using  phase  averages,  but  identical  results  can  be 
derived  using  equivalent  averages  over  phase  shift;  in  this  case 
the  averaging  is  one  over  a random  superposition  of  linear 
waves  and  is,  perhaps,  more  relevant  to  the  physical  situation). 
In  his  linear  analysis,  Landaihl  used  the  modified  action  equation 
accounting  for  weakly  nonconservative  effects,  as  given  in 
Eq.  (3^)*  The  actual  kinematics  are  therefore  describable  (to 
this  order)  by  the  "real"  dispersion  relation  JZ=ji.(K,  €X,  et), 
and  the  equation  for  wave  conservation  can  then,  within  the 
context  of  linear  theory,  be  used  to  predict  when  space/time 
singularities  in  the  wave  trajectories  can  be  expected  to  occur. 
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In  fact,  some  manipulation  shows  that  (Landahl,  1972)  focusing 
occurs  whenever 

(t)  ~^(V~dK  ^ 

is  satisfied,  where  Cq  is  the  critical  velocity  - 
In  the  steady  case,  the  critical  condition  reduces  to  the 
approach  toward  a zero  in  the  linear  group  velocity  (Landahl' s 
linear  theory  does  not  account  for  the  nonlinear  response  of 
the  supporting  medium).  This  condition  in  itself  implies 
increased  wave  amplitudes  on  account  of  "wave  tube"  convergence. 

If  this  slowing  down  is  additionally  accompanied  by  trapping, 
self-excited  waves  (as  determined  from  a local  application  of 
conventional  normal  mode  theory)  would  then  surely  amplify 
forever,  at  least  on  the  basis  of  linear  theory  (details  are 
given  in  Landahl,  1972).  One  would  then  surmise  that  the  sudden 
onset  of  this  kind  of  strong  instability  would  precipitate  a 
highly  nonlinear,  unsteady  readjustment  of  the  mean  flow,  in 
much  the  same  way  that  wave  trapping  and  amplification  in 
transonic  nozzles  lead  to  shock  formation  and  choking.  Of  course, 
Landahl* s analysis  is  a linear  one  and  therefore  does  not  account 
for  the  nonlinear  response  of  the  mean  flow.  But  the  physics  of 
the  present  study  is  in  the  same  spirit  of  Kantrowitz's  classic 
analysis,  in  which  the  small-perturbation  isentropic  equations 
are  used  to  describe  the  onset  of  the  instability  mechanism. 
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only  because  the  general  unsteady,  nonlinear  problem  is 
intractable.  Then,  as  in  Kantrowitz  (194?),  if  the  model  in 
fact  describes  the  true  hypothesized  instability  process,  one 
can  therefore  proceed  to  study  the  properties  of  the  shocks 
that  are  formed  in  the  latter  stages,  (Of  course,  one  can  study 
shocks  for  their  ovm  sake,  but  it  is  often  helpful  to  envision 
some  physical  generating  mechanism.  Some  of  the  ideas  proposed 
here  might  be  tested  by  numerical  integration  of  the  complete 
system  of  unsteady,  nonlinear,  high-order  modulation  equations 
subject  to  appropriate  initial  conditions.  Extreme  care  in  the 
analysis  must  be  taken,  however,  since  the  results  depend 
sensitively  on  the  delicate  balance  between  high-order  wave 
dispersion  and  diffusion,  not  to  mention  the  artificial  dispersive 
and  diffusive  effects  of  truncation  errors.) 

The  idea  of  a hypothetical  instability  mechanism  possibly 
leading  to  a final  configuration  with  discontinuous  end  states 
is  not  new,  and  some  inferences  can  at  least  be  drawn  from 
hydraulic  jump  formation  and,  in  vortex  breakdown,  they  are  at 
least  suggested  by  the  "finite  transition"  or  "conjugate  flow" 
theories  of  Benjamin  ( l%2  ),  This  latter  analysis  was  first 
to  demonstrate  the  possibility  of  an  infinitesimal  transition 
from  a supercritical  state  which  cannot  support  standing  waves, 
to  a subcritical  state  which  can  (but  the  theory  does  not  offer 
an  explanation  as  to  the  origin  of  the  discontinuity)!  the 
theory  essentially  postulates  in  an  ad  hoc  manner  the  physical 
necessity  for  the  existence  of  waves  on  the  subcritical  side,  so 
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as  to  maintain  the  continuity  of  total  momentum  flux.  We  can 

improve  on  Benjamin's  theory  by  considering  final  end  states 

that  are  essentially  inviscid  and  essentially  "organized"  in  the 

sense  that  a wave  description  is  still  applicable!  the  end 

state  of  the  "downstream"  portion  v/ould  consist  of  a highly 

enhanced  wave  riding  on  a mean  flow  of  reduced  energy,  formed 

from  an  "upstream"  portion  whose  v/ave  is  initially  infinitesimal 

in  amplitude  and  which  rides  on  a more  energetic  mean  flow.  We 

can  modify  Benjamin's  conjugate  flow  approach  (which,  in  effect, 

considers  only  the  mean  flow  in  developing  the  discontinuous 

solutions)  by  connecting  the  end  states  through  the  introduction 

of  two  more  free  parameters,  a wavenumber  K and  a wave  energy 

density  E of  a superimposed  wave  system  (however  small),  in 

addition  to  the  man  velocity  and  the  mean  cross-sectional 

distance.  In  the  context  of  hydraulic  jumps  (which  are  here 

treated  for  simplicity),  the  proposed  model  would  connect 

E,  K,  the  mean  height  h and  the  mean  speed  U by  conserving  across 

the  discontinuity  total  energy,  total  momentum,  total  mass  and 

waves  (i,e.,  we  fix  the  frequency,  for  a total  of  four  jump 

conditions).  This  model  thus  augments  Benjamin's  by 

systematically  accounting  for  the  presence  of  a wave  however 

small,  and  through  the  nonlinear  coupling,  establishes  whether 

or  not  conjugate  solutions  do,  in  fact,  exist  for  all  parameter 
A 

ranges.  Let  U be  the  depth-independent  "mass  transport 
velocity",  as  given  in  Phillips  (1969)»  which  already  accounts 
for  the  mass  flow  induced  by  the  wave.  For  simplicity,  we 
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examine  the  shallow  water  case,  so  that  the  wavenumber  K does 
not  explicitly  appear  in  the  conservation  laws  for  total  mass, 
momentum  or  energy.  The  results  v/e  are  to  derive,  of  course, 
are  valid  in  quasi-steady  problems  in  which  the  local  coordinates 
move  with  the  shock.  In  these  steady  coordinates,  the  governing 
equations  for  total  mass,  momentum  and  energy  integrate  to,  for 
a flat  bottom, 

l)U  ^ 7^ 

ifhu^  + JJh  E 


where  m,  P and  Q are  the  prescribed  (total)  mass,  momentum 

and  energy  fluxes,  P is  the  fluid  density,  h is  the  mean 

/ A 

depth  and  u is  the  rela-iive  velocity  U-Ug  (Ug  being  the  shock 
speed).  Let  us  introduce  the  nondimens ional  variables 

I = h/(^f^  ~ WK  , u ‘ ’ Vk 

so  that  the  above  conservation  laws  become,  respectively. 


hu.  f 


(104) 
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T 

flu} +e  = ~ ^ o (105) 

lju^+‘fuk^+fS2li''^£  -i-jEu^ ^3  - a > O 

If  E = 0 (bars  have  been  dropped  for  convenience)  and  Eq,  (106) 
is  deleted,  one  obtains  the  classical  hydraulic  jump  formulation 
in  which  mean  flow  mass  and  momentum  fluxes  are  fixed 
(Lamb,  1932).  However,  let  us  eliminate  u from  Eq.  (105) 
using  Eq.  (104)  to  produce  an  equation  for  E (in  terms  of  h 
and  P)  which,  in  turn,  is  used  in  Eq.  (106)  to  give  an  equation 
for  h.  Setting  H = h^  , a ninth  order  equation  is  obtained  as 

=o  (107) 

where  P,  Q and  H are  positive  and  all  roots  obtained  must  be 
such  that  E ^ 0 to  qualify  as  physically  realistic  conjugate 
solutions.  A rough  idea  for  the  type  of  jumps  typical  of 
Eqs.  (104-106)  is  obtained  by  comparison  with  the  results  of 
the  classical  formulation.  We  denote  by  "1"  and  "2"  the 
upstream  and  downstream  states  respectively,  and  require 
hg  > hj^  in  the  usual  way.  For  the  case  P = 2.12  (corresponding 
to  a Froude  number  of  3^  ) we  have  ^1  " ^*55  and 

h2  “ 1.10.  The  mean  flow  energy  flux  (in  the  classical  case) 
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Qjjj  = hu^  4-  4uh^  = 4h  + h“^  is  5.50  upstream  and  5.23  downstream, 
the  deficit  ratio  being  5.23/5.50  = 95.0^.  The  wave  theory  we 
propose  gives  similar  results,  but  here  the  loss  in  mean  flow 
energy  is  made  up  for  by  the  presence  of  radiating  waves.  For 
comparison’s  sake,  set  P = 2,12  and  Q = 5.50,  The  conjugate 
depths  are  now  obtained  as  h^  = 0,55  and  h2  = 1.03,  and 
calculations  show  = 0.0  while  ££  = 0,084,  In  other  words,  the 
waves  are  practically  nonexistent  on  Side  "1”,  but  on  Side  "2" 
they  acquire  a sudden  visibility.  The  ratio  of  downstream  to 
upstream  mean  energy  flux  is  92.1%  (as  opposed  to  95.0%  before), 
and  the  rest  goes  into  waves.  The  phase  velocity  can  be  shown 
to  decrease  across  the  jump  so  that,  since  the  frequency  is 
fixed,  the  wavenumber  increases.  The  sudden  visibility  of  the 
wave  is  enhanced  by  both  increased  wave  amplitude  and  wavenumber; 
the  classical  mean  flow  results  as  described  in  Lamb  (1932) 
are,  in  fact,  qualitatively  reproducible  by  the  above  wave 
model  without  special  appeal  to  turbulent  dissipation. 

More  general  results  are  obtained  by  examining  positive  E 
solutions  of  Eq,  (107).  There  is  a simple  way  to  construct  the 
entire  family  of  solutions,  and  this  is  seen  by  rewriting 
Eq,  (107)  in  the  form 


P 


(108) 


Thus,  every  value  of  the  positive  parameter  H defines  a 
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straight  line  in  the  P-Q  plane.  The  intersection  of  two  lines 
(or  more)  at  a point  (Pq'^o^  shows  that  there  exists  two  (or 
more)  conjugate  values  of  H at  that  point.  This  interpretation 
requires  us  only  to  study  the  slope  and  P- intercept  functions 
of  Eq,  (108) j not  every  point  on  a straight  line  in  the  P-Q 
plane  is  physically  significant,  though,  since  it  is  necessary 
for  the  corresponding  wave  energy  to  be  positive.  The  simplest 
criterion  for  this  is  found  from  Eq.  (105),  that  is, 

^ ^ V 'p-  ^ O 

Thus,  for  E ^ 0,  it  is  necessary  that  P ^ P*  = -t-  H“^.  A 

useful  presentation  of  some  results  is  shown  in  Figure  1. 
Various  sets  of  are  considered  such  that  h2/hj^  = 2.25, 

so  that  we  in  fact  assume  (H^,H2)  to  be  a conjugate  pair. 

Then,  each  H determines  a straight  line  in  the  P-Q  plane,  with 
the  intersection  point  Pq»Qq  determined  from  the  solution  of 
two  simultaneous  linear  equations  derived  from  Eq.  (108). 
Repeating  the  process  for  various  h's  generates  different 
solutions,  as  shown  in  Figure  1 (different  values  of  h^/h-^  were 
also  considered,  and  the  same  qualitative  results  were 
obtained).  The  results  show  that  every  point  along  the  solid 
portion  of  the  line  determines  two  conjugate  solutions i the 
jump  in  h from  h^^  to  h^  is  always  accompanied  by  an  increase  in 
the  wave  energy  density  and  a decrease  in  the  mean  flow  energy 
flux  Qju  previously  defined,  again  agreeing  qualitatively  with 
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classical  theory,  but  without  appeal  to  turbulent  dissipation. 

Each  point  along  the  solid  curve  here  possesses  two  roots  H 
with  positive  E’s,  while  the  dashed  portion  indicates  the 
existence  of  one  negative  root}  in  this  latter  case  the  solution 
is  unique  and  conjugate  solutions  are  not  possible.  It  follows 
that  conjugate  solutions,  at  least  in  the  case  of  shallow 
water  gravity  waves,  are  possible  only  for  sufficiently  large 
P's  or  Q’ s. 

A second  kind  of  discontinuous  solution  is  typified  by  the 

flow  of  gravity  waves  over  infinitely  deep  constant  speed 

currents.  Because  the  current  in  this  overly  simplified  model 

does  not  vary  with  depth,  jumps  in  the  mean  speed  must  be 

disallowed,  for  they  would  otherwise  imply  infinite  changes  in 

mean  energy.  The  wave  in  this  sense  decouples  from  the  mean 

flow  and  the  only  discontinuities  that  form  are  those  in  the 

wave  parameters  themselves  (jumps  in  for  y-dependent  mean 

A 

flows  can  conceivably  be  analyzed  in  terras  of  jumps  in  U as 
given  in  Eq,  (51)i  but  it  is  unclear  as  to  whether  the  corresponding 
shape  function  remains  invariant  across  the  discontinuity).  By 
varying  the  parameter  a family  of  related  conjugate  solutions 
can  be  obtained  and  these  solutions  will  form  the  present  focus 
of  attention.  It  will  be  helpful  here  to  envision  surface  or 
plasma  waves  in  steady  flow  and,  for  reasons  previously  indicated, 
discontinuities  across  which  frequency  ( (/Jq  ) momentum 
flux  ( ) are  fixed  (the  modifications  needed  to  consider  wave 
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breaking  are  not  considered  here).  Some  analysis  shows  (this 
will  be  verified  later)  that  when  the  (real)  frequency  is 
expanded  in  the  wave  action  density  A,  that  is, 

V-„K*  F(K)  +AGCK)  = 

the  corresponding  expression  for  momentum  flux  becomes 

Solutions  can  be  readily  obtained  by  eliminating  between 
these  two  equations,  tabulating  A as  the  wavenumber  K is  varied, 
then  solving  for  U_  from  either  of  the  above  equations,  and 
finally  solving  for  the  momentum  velocity  defined  by 

^ K y 

The  resulting  equations  are  valid  for  both  waves  advancing  and 
receding  with  respect  to  the  mean  flow,  of  course,  provided 
appropriate  sign  changes  in  P(K)  and  G(K)  are  made.  Only  real 
solutions  of  A with  k ^ 0 are  admissible,  and  hence,  multivalued 
solutions  may  or  may  not  exist  for  all  values  of  the  parameter 
For  illustrative  purposes  let  us  consider  two  weakly 
nonlinear  irrotational  examples,  first  gravity  waves  and  second, 
capillary  waves.  In  the  case  of  gravity  waves,  the  appropriate 
dispersion  functions  F ■ g^K^  and  G = where,  again,  g is  the 

acceleration  due  to  gravity.  The  barred  nondimensional  variables 
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4^ 

^ yyj  ^ 

U = jM  U„  c = f V./  f - 2S&  - jj- 

are  introduced,  and  we  shall  assume  without  loss  of  generality, 

that  OU  >0  throughout.  Figure  2 shows  the  complete  family  of 
o 

p,C  ^ 0 "positive  moving”  waves  consisting  of  both  waves 
advancing  and  receding  with  respect  to  the  mean  flow,  while 
Figure  3 shows  the  complete  solution  for  "negative  moving” 
waves  with  p,C  < 0,  For  instance,  from  Figure  2,  only  one 
wave  solution  exists  when  U<i<  -1  but  three  solutions  exist  for 
U ^ 0,  say,  corresponding  to  a given  flux  p.  Figure  3 shows  how 
for  positive  U,  two  wave  solutions  can  be  found  for  a given 
negative  flux  p < 0 ( corresponding  to  a "strong"  high-amplitude 
and  a "weak"  low-amplitude  solution  across  which  energy  flux 
decreases),  but  solutions  for  sufficiently  negative  U’s  cannot 
be  found  (this  may  point  to  the  importance  of  unsteady  effects). 
The  "singular  point"  occurring  at  U = -1  approximately 
corresponds  to  the  position  determined  from  linear  theory  as 
that  producing  a zero  in  the  linear  group  velocity.  For 
capillary  waves  (where  F * and  G = -K^/8)  the  (barred) 

nondiaensional  variables 
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with  surface  tension  T and  nondimensional  aomentum  flux 

o 

3 = m /T  can  be  defined.  Let  us  here  assuae  UJ  0 
0 0 o 

throughout.  For  various  specified  values  of  positive  s,  the 
complete  solutions  (consisting  of  both  waves  advancing  and 
receding  with  respect  to  U ) are  plotted  in  Figure  4;  negative 
moving  solutions  with  s ^ 0 are  shown  in  Figure  5.  From 
Figure  4 for  positive  moving  waves,  multivalued  solutions 
corresponding  to  various  fixed  values  of  s exist  for  sufficiently 
negative  "U's,  but  for  positive  U’s  the  solutions  are  single- 
valued; similar  considerations  hold  for  negative  moving  waves, 
as  is  evident  from  Figure  5.  Again,  the  U position  about  which 
the  multi-valued/single-valued  character  of  the  problem  pivots 
is  determined  approximately  by  the  value  of  U for  which  the 
linear  group  velocity  is  zero.  On  Figures  2-5  the  nondimensional 
wavenumber  distributions  are  also  shown. 

In  practice  the  jumps  discussed  above  are  formed  on 
currents  that  vary  slowly  in  space  and  time,  and  it  is 
important  here  to  note  that  the  above  results  are  applicable 
as  well  in  quasi- stationary  coordinates  moving  with  the  shock. 

One  question  of  fundamental  importance  concerns  the  observability 
of  the  hypothesized  shocks,  and  this  can  be  answered  in  part  by 
examining  their  stability  with  respect  to  self-induced 
disturbances.  In  the  simple  analysis  that  follows,  high-order 
modulation  terms  (that  in  reality  allow  or  disallow  the  shock 
formation)  will  be  ignored;  we  in  effect  assume  the  existence  of 
a shock  and  study  its  stability  as  would  be  affected  by  low-order 
inhomogeneities  in  the  medium.  Two  types  of  shocks  can  i''.:  • 
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general  treatment,  and  they  have  already  been  introduced 
in  the  above  analyses.  Type  A shocks,  like  those  just  treated, 
conserve  wave  momentum  and  waves,  while  the  wave  energy  flux 
decreases  across  the  discontinuity.  The  mean  flow  here  enters 
only  parametrically  and  it  formally  uncouples  from  the  wave 
system.  Type  B shocks  are  defined  as  those  in  which  the  mean 
flow  responds  dynamically  to  changes  in  the  wave  parameters, 
for  example,  as  in  the  generalized  hydraulic  jumps  considered 
previously  for  water  of  finite  depth i these  by  definition  invoke 
the  constancy  of  total  mass,  total  momentum,  total  energy  and 
frequency. 

We  shall  consider  first  Type  A shocks.  It  is  convenient 
to  assume  a perturbation  Lagrangian  ^ in  the  form  <l^=AJl-H(K,A) 
as  viewed  from  a coordinate  system  moving  with  the  mean  flow 
U_  I henceforth  all  subscript  m's  will  be  deleted,  (Previously 

Jn 

wo  assumed  an  expanded  H function  written  in  powers  of  A,  with 
H = AF(K)  + iA^G(K),  The  frequency  following  the  flow  U was 
thus  determined  from  ■ 0 as  F+  AG,  implying  a 

ground-fixed  frequency  of  the  form  = UK+  F+-AG,  The 

corresponding  momentum  flux  was  determined  from  UKA  + X -kX^,  ) 
To  begin  the  analysis,  we  note  that  an  unsteady  shock  with 
speed  Ug  • satisfies  the  conditions  for  frequency  and 

momentum  flux  continuity  in  the  following  respective  forms 


Cv-v,) A,K “ ••• 


(no) 
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where  the  *1"  and  "2"  subscripts  refer  to  upstream  and 
downstream  portions  of  the  discontinuous  solution.  Let  us 
consider  the  basic  problem  in  which  a stationary  shock  is 
displaced  a small  amount  \ , A,  K and  U can  thus  be  expanded 
by  Taylor  series  about  the  steady  solution  corresponding 
to  I - 0,  noting  that  the  shock  speed  Ug(T)  = f(T)  is  a 
quantity  of  0(  ^ )•  For  example,  the  suggested  expansions  give 

A ^ ^0*^0  ^ function  S(A),  S(A)  = ^ ^ ^A^'^o 

To  0(1)  we  have  the  following  frequency  and  momentum  flux  jump 

conditions 


(111) 

(112) 

the  weakly  nonlinear  forms  of  which  were  used  in  the  previous 
two  examples,  "o"  subscripts  here  indicating  the  basic  mean 
state.  The  0(  ^ ) equations  lead  to  the  equations 
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1 

I 


(11^) 


where  primes  denote  1 derivatives.  As  in  Kantrowitz  (19^7) i a 
number  of  growth  rates  exist,  each  of  which  corresponds  to 
various  permissible  values  of  the  shock  speed  U^,  and  are  to 
be  separately  investigated.  Equations  (113)  and  (114)  can  be 
simplified  by  invoking  wave  and  action  conservation  on  either 
side  of  the  discontinuity  (but  not  across  it).  On  this  basis 
one  finds  that  each  of  the  bracketed  quantities  in  Eq.  (113) 
separately  vanishes,  implying  a neutrally  stable  shock  solution. 
The  second  solution,  obtained  from  Eq,  (114),  is  remarkably 
simple I 


I » Vj^o  I (115) 

Now  the  mean  speed  U(X,T)  in  the  present  model  is  prescribed, 

and  to  fulfill  the  slowly  varying  assumption,  the  derivative 

U1  = Ui  s (with  respect  to  X in  a steady  coordinate 

2,0  1,0  (0) 
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system)  must  be  continuous.  Thus,  Eq.  (115)  deduces  to 

• 

^ 3 £ , exactly  the  same  condition  for  shock  stability 

transonic  nozzles!  A shock  is  therefore  stable  or  unstable 
accordingly  as  ® ^ simple  criterion 

showing  how  the  shock  stability  is  wave  amplitude  and  wavenumber 

independent.  Results  for  U*  ^ 5^  U!  can  also  be  discussed,  of 

2,0  1,0 

course,  but  this  would  appear  to  be  of  limited  usefulness 
since  the  slowly  varying  assumption  would  require  continuity 
of  second  derivatives  of  U as  well.  In  a uniform  medium,  the 
shocks  considered  here  would  be  neutrally  stable. 

Type  B shocks  (which  allow  for  the  dynamic  response  of  the 
mean  parameters)  can  be  treated  by  considering  in  addition  to 
iktUJ  ,Yi)  the  triad  ).  In  water  waves,  for  example, 

P and  ^ would  be  the  mean  height  and  mean  speed,  respectively, 
and  »P*A,K)  would  be  a Bemoulll-type  "dispersion 

relation"  for  the  pressure.  We  take  the  Lagrangian  in  the  form 


* A(jU  F,^) 


(116) 


following  Hayes  (1973)i  where  A = A (X,T)  describes  the 
iiUiomogeneity,  The  total  energy  of  the  system  is  H 
and  ^ ■ 0 defines  the  dispersion  relation  Cx/=  H^.  A second 
dispersion  relation  arises  from  taking  - 0,  that  is, 

JT  * Hp  , Thus  there  remains  foxxr  unknowns,  for  which  we  invoke 
the  conservation  laws  for  waves,  total  mass,  total  energy  and 
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total  monentuR.  They  are,  respectively, 

^ ?uj 

at  ax  * 0 

(117a) 

(118a) 

(119a) 

(120a) 

In  a coordinate  system  moving  with  the  shock, 

jump  conditions  become 

a 

the  respective 

-ir^K  * Car,st«H-t 

(117b) 

i 

(118b) 

(119b) 

(120b) 

As  beforst  a basic  steady  discontinuous  solution  is  assuaed, 
and  the  effect  of  sllfhtly  displacing  the  discontinuity  is 
exaained.  An  analysis  siailar  to  that  Just  coapleted  for 
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Type  A shocks  leads  to,  for  the  0(1)  heirarchy  of  equations, 
the  juap.  conditions  for  the  basic  steady  flowi 


while  the  0(  ^ ) equations  are,  respectively, 

j ( X -(  , 

c ; “ + HmA' 

{ - c x-c ),  , 


(117c) 

(118c) 

(119c) 

(120c) 


(117d) 


(ll8d) 
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l -ivlv  ‘ 

+ H^CHjicIc'  +tlf^A'  + Wff('  ■<-WfrT'-^Wf;,X) 

■''H^^Hpi^fC  ■'■HmA  +H^pT +Hf;iA,^ 

C - ^’-  I (120d) 

'( (ftjoA 

( )”  A(H^K'+H^A'  + ri*^^'+H/(fT 

*?(Hn<'-^  Hp^^'+  HT/'-*-HTr7’'’-M0 
+^r(HKKr +«»,  A' Hk-^  A h„p'+Hic(  v; 

■'*^f(l^fK^' 

-H/ 

The  bracketed  quantities  in  Eqs.  (117d)(  (118d)t  (119d)  and 
(120d)(  of  course,  can  be  siaplified  by  applying  the  steady 
differential  equations  (for  example,  the  wave  action  equation 
dil^dx  = 0,  the  wave  conservation  law  PH^/^x  = 0,  PH^/px  a 0 
for  aass  conservation,  and  the  consistency  equation 
P^/Pt  ♦ ^y/dx  a 0 or  3Hp/a  X = 0 ).  On  substitution,  the 
amplification  rates  in  the  first  three  equations  can  be  shown  to 
be  identically  zero,  implying  neutral  stability.  However,  the 
result  for  Eq.  (120d)  is  nontrivial  and  leads  to 


101 


(121) 


Thus,  if  the  inhomogeneity  2,'  ''banishes,  the  shock  is  neutrally 
stable.  Some  applications  of  this  formula  cam  be  found  in 
Chin  (1976),  but  the  qualitative  results  are  similar  to  those 
obtained  for  Type  A shocks. 

The  basic  ideas  behind  the  instability  theory  advanced 
in  this  section  can  be  briefly  summariaedi  an  infinitesimally 
small  disturbance  wave  riding  on  a slowly  varying  mean  flow,  on 
account  of  space-time  inhomogeneities,  becomes  trapped  in  a 
region  that  on  the  basis  of  local  normal  mode  theory  is  unstable. 
The  rapid  amplification  of  the  trapped  wave  then  precipitates 
a sudden  nonlinear,  unsteady  readjustment  of  the  mean  flow  that 
is  locally  viewed  as  a "generalized  hydraulic  jump"  involving 
both  mean  flow  and  wave  parameters.  The  actual  observability 
of  the  hypothesized  "shock"  mechanism  would,  of  course,  depend 
to  a great  extent  on  the  all-important  effect  of  high-order 
dispersive  and  diffusive  terms  as  well  as  on  the  effect  of 
low-order  shock  stability.  These  ideas  on  the  origin  of 
"catastrophic  transitions"  are  at  the  present  time  speculative 
in  nature,  but  they  do  appear  to  be  motivated  by  plausible 
physical  arguments  and  by  gas-dynamical  analogies  1 they 
explain  the  formation  of  the  discontinuities  assumed  in,  for 


t 


example,  the  classical  treatment  of  hydraulic  jumps  and 
Benjamin's  conjugate  flow  theory  for  vortex  breakdown.  Because 
the  nonlinearities  of  a physical  system  are  particularly 
important  in  its  detailed  description,  the  possibility  of 
deducing  general  results  seems  remote.  Application  of  the 
generalized  high-order  Whitham  theory  developed  in  this 
paper,  for  example,  to  hydraulic  jump  foraiation,  vortex 
breakdown  or  boundary  layer  transition  in  this  regard  would  be 
extremely  interesting.  One  immediate  need,  perhaps,  is  a 
nonlinear  wave  trapping  theory  paralleling  Landahl's  for 
linear  flows i some  initial  work  has  been  carried  out  in 
Chin  (1976),  but  much  more  work  needs  to  be  done. 

V.  Summary. 

A theory  has  been  constructed  for  the  description  of 
wavetrains  slowly  modulated  in  apace  and  time  and  which  includes 
the  effects  of  high-order  diffusive  and  dispersive  corrections, 
modal  terms,  nonconservative  terms,  and  amplitude  and  frequency 
dispersion.  A number  of  examples  have  been  considered  leading 
to  many  formulas  useful  in  direct  applications  of  the  theory, 
and  various  new  ideas  in  stability  have  been  discussed. 
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FIGURE  2 nonlinear  gravity  waves,  CU  > 0 AND  C>0. 
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FIGURE  3A 


figure  4 

Honltnear  Capillary  Waves,  0 "^>0. 

Wavenumber  Distribution  Plotted  on  Curve. 
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FIGURE  5 

Nonlinear  Capillary  Waves,  0 and  C<(  0. 

Wavenumber  distribution  plotted  on  curve. 
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